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FOREWORD 

The  purpose  of  this  study  is  to  survey  the  analytical  linear  solutions  for 
a  thin  cylindrical  shell  loaded  by  point  forces  and  compare  them  to  various 
finite  element  models.  The  shell  theories  used  to  quantify  maximum  central 
displacement  are  the  so-called  "technical"  or  "engineering"  theories.  A 
computer  program  was  developed  to  automate  the  procedure,  and  a  variety  of  plots 
for  various  radius-to-thickness  and  length-to-radius  ratios  vas  obtained.  The 
finite  element  solutions  were  computed  using  general  purpose  computer  programs 
STAGS  and  ABAQUS.  The  agreement  between  the  analytical  and  the  finite  element 
results  is  quite  good. 

This  work  was  sponsored  by  NAVSEA  (63R)  through  the  Undersea  Weaponry 
Warheads,  and  Fuzes  Block  Program. 

Approved  by: 

KURT  F.  MUELLER,  Head 

Energetic  Materials  Division 


»■  a*  -ar  ,"w»..-bt  ,  M  rji  "ai 


NSWC  TR  84-342 


CONTENTS 

Chapter  Page 


1  INTRODUCTION .  I 

2  GOVERNING  EQUATIONS  .  3 

3  ANALYTICAL  SOLUTIONS .  5 

ODQVIST  SOLUTION .  5 

BIJLAARD'S  SOLUTION  .  6 

MIZOGUCHI'S  SOLUTION .  6 

YUAN-TING’S  SOLUTION .  7 

4  NUMERICAL  SOLUTIONS  .  11 

POINT  LOAD  CASE .  11 

PINCHED  LOAD  CASE .  12 

5  COMPARISON  OF  NUMERICAL  AND  ANALYTICAL  SOLUTIONS .  13 

POINT  LOAD  CASE .  13 

PINCHED  LOAD  CASE .  14 

CONCLUSION .  13 

REFERENCES .  49 

NOMENCLATURE .  55 

Appendix  Page 


A  SOLUTION  METHODS  .  A-l 

B  BIJLAARD'S  SOLUTION  BY  DOUBLE  FOURIER  SERIES .  B-l 

C  SYMMETRY  OF  VLASOV  SHELL  THEORY  .  C-l 


■  j  w  r  i 


Figure 


NSWC  TR  84-342 


ILLUSTRATIONS 


CYLINDRICAL  SHELL  SIMPLY  SUPPORTED  AT  BOTH  ENDS  AND  LOADED  BY  A 
RADIALLY  INWARDS  CONCENTRATED  FORCE  P.  SHELL  IS  OF  RADIUS  R, 
THICKNESS  h,  AND  LENGTH  L  . 


PINCHED  CYLINDRICAL  SHELL  BY  INWARDS  CONCENTRATED  FORCES  P, 
SHELL  IS  OF  RADIUS  R,  THICKNESS  h,  AND  LENGTH  L  . 


NONDIMENS IONAL  PARAMETER  EHW/P  FOR  A  SIMPLY  SUPPORTED  CYLINDRICAL 
SHELL  SUBJECT  TO  A  RADIALLY  INWARDS  CENTRALLY  LOCATED 
CONCENTRATED  FORCE  P  (BY  ODQVIST  SOLUTION)  PLOTTED  VERSUS  R/L 
FOR  VARIOUS  R/H  RATIOS . 

NONDIMENS IONAL  PARAMETER  EHW/P  FOR  A  SIMPLY  SUPPORTED  CYLINDRICAL 
SHELL  SUBJECT  TO  A  RADIALLY  INWARDS  CENTRALLY  LOCATED 
CONCENTRATED  FORCE  P  (BY  BIJLAARD  SOLUTION)  PLOTTED  VERSUS  R/L 
FOR  VARIOUS  R/H  RATIOS . 

NONDIMENS IONAL  PARAMETER  EHW/P  FOR  A  SIMPLY  SUPPORTED  CYLINDRICAL 
SHELL  SUBJECT  TO  A  RADIALLY  INWARDS  CENTRALLY  LOCATED 
CONCENTRATED  FORCE  P  (BY  MIZOGUCHI  SOLUTION)  PLOTTED  VERSUS  R/L 
FOR  VARIOUS  R/H  RATIOS . 

NONDIMENS IONAL  PARAMETER  EHW/P  FOR  A  SIMPLY  SUPPORTED  CYLINDRICAL 
SHELL  SUBJECT  TO  A  RADIALLY  INWARDS  CENTRALLY  LOCATED 
CONCENTRATED  FORCE  P  (BY  YUAN-TING  SOLUTION)  PLOTTED  VERSUS  R/L 
FOR  VARIOUS  R/H  RATIOS . 

NONDIME NS IONAL  PARAMETER  EHW/P  FOR  A  RADIALLY  INWARDS  PINCHED 
CYLINDER  WITH  FREE  ENDS  (BY  INEXTENS IONAL  SOLUTION)  PLOTTED 
VERSUS  L/R  FOR  VARIOUS  R/H  RATIOS  . 


NONDIMENS IONAL  PARAMETER  EHW/P  FOR  A  SIMPLY  SUPPORTED  CYLINDRICAL 
SHELL  SUBJECT  TO  A  RADIALLY  INWARDS  CENTRALLY  LOCATED 
CONCENTRATED  FORCE  P  (BY  BIJLAARD  SOLUTION)  PLOTTED  VERSUS  R/L 
FOR  VARIOUS  R/H  RATIOS . 

NONDIME NS IONAL  PARAMETER  EHW/P  FOR  A  SIMPLY  SUPPORTED  CYLINDRICAL 
SHELL  SUBJECT  TO  A  RADIALLY  INWARDS  CENTRALLY  LOCATED 
CONCENTRATED  FORCE  P  (BY  MIZOGUCHI  SOLUTION)  PLOTTED  VERSUS 
R/L  FOR  VARIOUS  R/H  RATIOS . 


ft 


$ 


u 


24 


/wjiidbnity  Codes 

Avail  aivd/or 
Special 


/?-/ _ 


NSWC  TR  84-342 


ILLUSTRATIONS  (Cont.) 


F igure  Page 

10  NONDIMENS ION AL  PARAMETER  EHW/P  FOR  A  SIMPLY  SUPPORTED  CYLINDRICAL 

SHELL  SUBJECT  TO  A  RADIALLY  INWARDS  CENTRALLY  LOCATED 

CONCENTRATED  FORCE  P  (BY  YUAN-TING  SOLUTION)  PLOTTED  VERSUS 

R/L  FOR  VARIOUS  R/H  RATIOS .  25 

11  NONDIMEN SIONAL  PARAMETER  EHW/P  FOR  A  RADIALLY  INWARDS  PINCHED 

CYLINDER  WITH  FREE  ENDS  (BY  INEXTENS IONAL  SOLUTION)  PLOTTED 

VERSUS  L/R  FOR  VARIOUS  R/H  RATIOS .  26 

12  SUPERIMPOSED  PLOTS  OF  ANALYTICAL  SOLUTIONS  BY  FOUR  METHODS  AND 

FOUR  FINITE  ELEMENT  DISCRETIZATIONS  OF  EHW/P  VERSUS  L/R  FOR 

R/H  =  100 .  27 

13  SUPERIMPOSED  PLOTS  FOR  A  PINCHED  CYLINDER  BY  TIMOSHENKO'S 

INEXTENS IONAL  SOLUTION  AND  TWO  FINITE  ELEMENT  DISCRETIZATIONS  .  28 


M 


14  SIMPLY  SUPPORTED  CYLINDER  (9  *  180°)  SUBJECTED  TO  RADIAL 

CONCENTRATED  LOAD  (BY  STAGS) .  29 

15  SIMPLY  SUPPORTED  CYLINDER  U  -  72",  R  -  100",  9  =90°)  SUBJECTED 

TO  RADIAL  CONCENTRATED  LOAD  VIEWED  FROM  (100",  100",  100")  (BY 
ABAQUS) .  30 


B-l  CYLINDRICAL  SHELL  SIMPLY  SUPPORTED  AT  BOTH  ENDS  AND  LOADED  BY  A 
RADIALLY  INWARDS  CONCENTRATED  FORCE  P.  SHELL  IS  OF  RADIUS  R, 
THICKNESS  h,  AND  LENGTH  L .  B”2 

B-2  SIMULATION  OF  SIMPLY  SUPPORTED  CONDITION  BY  CONSIDERING  AN 

INFINITE  NUMBER  OF  POINT  LOADS  P  AT  A  DISTANCE  L  APART .  B-8 


iv 


TABLES 


Table 

1  GOVERNING  DIFFERENTIAL  EQUATIONS  USED  FOR  ANALYTICAL 

SOLUTIONS  . 

2  SUMMARY  OF  SOME  PRACTICAL  METHODS  . 

3  GEOMETRICAL  AND  MATERIAL  PROPERTIES  OF  ANALYZED  MODELS  . 

4  RESULTS  BY  STAGSC  COMPUTER  PROGRAM  FOR  A  SIMPLY  SUPPORTED  SHELL 

(V«W=RU=*0)  SUBJECT  TO  A  CENTRALLY  LOCATED  INWARD  CONCENTRATED 
FORCE  . 

5  RESULTS  BY  STAGSC  COMPUTER  PROGRAM  FOR  A  SIMPLY  SUPPORTED  SHELL 

(V-W“RU=-0)  SUBJECT  TO  A  CENTRALLY  LOCATED  INWARD  CONCENTRATED 
FORCE  . 

6  RESULTS  BY  STAGSC  COMPUTER  PROGRAM  FOR  A  PINCHED  SHELL  (BOUNDARY 

CONDITIONS  AT  ENDS  RV-RW=*0) . 

7  CONVERGENCE  STUDY  BY  STAGSC  COMPUTER  PROGRAM  OF  MODEL  NO.  1 

WITH  CONTINUITY  END  CONDITIONS  . 

8  CONVERGENCE  STUDY  BY  STAGSC  COMPUTER  PROGRAM  OF  MODEL  NO.  1 

WITH  SIMPLY  SUPPORTED  ENDS  . 

9  RESULTS  BY  ABAQUS  COMPUTER  PROGRAM  FOR  A  SIMPLY  SUPPORTED  SHELL 

(V=W=*RU“0)  SUBJECT  TO  A  CENTRALLY  LOCATED  INWARD  CONCENTRATED 
FORCE  . 

10  RESULTS  BY  ABAQUS  COMPUTER  PROGRAM  FOR  A  SIMPLY  SUPPORTED  SHELL 

(V<*W-RU=*0)  SUBJECT  TO  A  CENTRALLY  LOCATED  INWARD  CONCENTRATED 
FORCE  . 

11  RESULTS  BY  ABAQUS  COMPUTER  PROGRAM  FOR  A  PINCHED  SHELL  (BOUNDARY 

CONDITIONS  AT  ENDS  RV-RW-O)  . 

12  PARAMETER  EWH/P  FOR  VARIOUS  R/H  AND  R/L  RATIOS  FOR  A  PINCHED 

CYLINDRICAL  SHELL  BY  TIMOSHENKO'S  INEXTENS IONAL  SOLUTION  .  .  .  . 

13  PARAMETER  EWH/P  FOR  VARIOUS  R/H  AND  R/L  RATIOS  FOR  A  SIMPLY 

SUPPORTED  CYLINDRICAL  SHELL  SUBJECT  TO  A  CENTRALLY  LOCATED 
INWARD  POINT  LOAD  BY  FOUR  ANALYTICAL  METHODS  . 


NSWC  TR  84-342 


TABLES  (Cont.) 

Table  Pa; 

14  PARAMETER  EWH/P  FOR  VARIOUS  R/H  AND  R/L  RATIOS  FOR  A  PINCHED 

CYLINDRICAL  SHELL  BY  TIMOSHENKO'S  INEXTENSIONAL  SOLUTION  ....  4< 

15  PARAMETER  EWH/P  FOR  VARIOUS  R/H  AND  R/L  RATIOS  FOR  A  SIMPLY 

SUPPORTED  CYLINDRICAL  SHELL  SUBJECT  TO  A  CENTRALLY  LOCATED 

INWARD  POINT  LOAD  BY  FOUR  ANALYTICAL  METHODS . 45 

16  PARAMETER  EWH/P  FOR  VARIOUS  R/H  AND  R/L  RATIOS  FOR  A  PINCHED 

CYLINDRICAL  SHELL  BY  TIMOSHENKO'S  INEXTENSIONAL  SOLUTION  ....  46 

17  PARAMETER  EWH/P  FOR  VARIOUS  R/H  AND  R/L  RATIOS  FOR  A  SIMPLY 

SUPPORTED  CYLINDRICAL  SHELL  SUBJECT  TO  A  CENTRALLY  LOCATED 
INWARD  POINT  LOAD  BY  FOUR  ANALYTICAL  METHODS 


i 

1 

\ 

A 

A 


NSWC  TR  84-342 


CHAPTER  1 
INTRODUCTION 


The  validation  of  finite  element  computer  programs  for  various  config¬ 
urations,  loads,  and  material  properties  is  of  great  interest  to  the  Navy,  which 
is  applying  these  methods  on  a  regular  basis. 

Structural  statics  problems,  depending  on  the  loading,  material  properties, 
and  geometry  changes  during  deformation  can  be  classified  as  geometrically 
linear  or  nonlinear,  elastic,  plastic  or  elastoplastic  with  small  or  large 
deflections  and/or  rotations.  For  convenience  we  will  refer  to  three  ranges  of 
interest,  i.e.,  (1)  the  small  deflection  elastic  material  regime,  (2)  the  large 
deflection  and/or  rotations,  fully  plastic  material  regime,  and  (3)  an 
intermediate  range  with  elastoplastic  properties  and  deformations  ranging  from 
small  to  large.  It  is  clear  that  all  three  ranges  must  be  validated. 

This  report  addresses  two  specific  structural  problems,  and  focuses  on  the 
first  and  simplest  range  of  interest,  i.e.,  the  small-deflection  linear  regime. 
The  two  problems  are: 

1.  A  circular  cylindrical  shell  simply  supported  at  both  its  ends  and 
subjected  to  a  centrally  located  inward  concentrated  static  load. 

2.  A  self-equilibrating  circular  cylindrical  shell  subjected  to  two  inward 
pinching  static  loads  at  equal  and  opposite  radial  positions. 

The  report  also  presents  a  comparison  of  finite  element  solutions  obtained 
using  both  the  STAGS^  and  ABAQUS^  computer  codes  with  analytical  solutions 
found  in  the  literature  for  thin  shells . 3 « 5 > 7  ,8-12  The  analytical  methods 
are  derived  from  the  so-called  "engineering"  or  "technical"  theories  of  thin 
shells,  which  differ  in  the  complexity  and  type  of  results  that  can  be  generated 
from  solutions  based  on  three-dimensional  elasticity  theory  methods.  The  two 
loading  cases  (1)  and  (2)  are  illustrated  in  Figures  1  and  2  and  labeled  as  the 
"point  load  case"  and  the  "pinched  load  case,"  respectively.  Other  loading  and 
boundary  conditions  are  examined  in  an  extensive  review  of  the  pertinent 
literature  (References  3  to  80),  some  of  which  is  included  as  Appendix  A. 

Finally  a  variety  of  plots  for  various  radius-to-thickness  and  length-to-radius 
ratios  were  obtained.  They  can  be  used  for  design  purposes  without  recourse  to 
a  computer,  since  they  are  given  for  nondimens ional  variables  and  a  range  of 
interest  to  marine  structures. 

Future  studies  will  address  the  same  two  problems  for  the  second  range  of 
interest,  i.e.,  for  a  fully  plastic  shell  treated  as  if  made  of  a  "rigid" 
plastic  material. 
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CHAPTER  2 

GOVERNING  EQUATIONS 


The  analytical  solutions  with  which  the  numerical  results  are  later 
compared,  follow  the  works  of  Odqvist,-^  Bijlaard,-^  Mizoguchi,-^ 

Yuan-Ting, ^  and  Timoshenko.^  These  authors  (except  Timoshenko  in  this 
case)  developed  partial  differential  equations  governing  the  response  of  the 
cylinder  under  the  loads  of  interest,  that  represent  simplifications  of  the  more 
general  detailed  governing  equations  as  represented  in  the  works  of  Fluegge.^ 

A  detailed  account  of  the  development  of  some  of  these  equations  employed  in 
this  report  is  presented  in  Appendix  B. 

Table  1  summarizes  the  derived  partial  differential  equations  used  in  this 
study,  and  on  which  the  computed  theoretical  solutions  are  based.  Table  2  is  an 
additional  summary  and  supplements  Table  1.  Timoshenko's  solution  was  developed 
by  employing  the  principle  of  virtual  work  and  hence  no  differential  equation 
was  directly  solved  and  given  in  Table  1.  Further  details  on  all  the  solutions 
can  be  found  in  Appendix  B  and  the  references  cited  therein. 
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CHAPTER  3 


ANALYTICAL  SOLUTIONS 


This  portion  of  the  report  deals  with  the  employed  analytical  solutions. 

In  Odqvist's,  Bijlaard's,  Mizoguchi's,  and  Yuan-Ting's  case  we  start  from  the 
equations  of  equilibrium  for  an  infinitesimal  shell  element,  introduce  the 
strain-displacement  equations,  and  follow  the  approach  outlined  in  Appendix  B  to 
reduce  the  equations,  to  a  partial  differential  equation  in  one  variable  after 
neglecting  certain  terms. 

Here  it  is  sufficient  to  say  that,  depending  on  the  type  of  assumptions 
made  with  respect  to  the  thickness-to-radius  ratio  (h/R) ,  the  resulting 
equations  can  be  simplified  accordingly.  Consequently,  we  can  have  on  the  one 
hand  the  full  Fluegge3  equations,  which  are  the  most  complicated,  while  on  the 
other  hand  the  Donnell3^*  equations,  which  are  the  most  simplified  consistent 
version.  We  also  note  that  the  even  simpler  equation  given  by  Schorer3^  and 
Odqvist30  will  be  used  for  numerical  comparisons  later  on.  In  addition  there 
are  other  variations  with  respect  to  complexity.  Further  discussion  can  be 
found  in  Appendices  A,  B,  and  C  and  the  references  therein. 

Figure  1  displays  the  point  load  case.  The  cylinder  is  loaded  at  midlength 
(L/2)  by  an  inwards  concentrated  force  P  and  is  simply  supported  at  its  two 
extreme  ends  x  =  0  and  x  *  L. 

Figure  2  displays  the  pinched  load  case.  The  cylinder  is  loaded  at 
midlength  (L/2)  by  two  equal  and  opposite  inwards  directed  forces  P  at  the 
angular  locations  q>  *  0  and  tp  =  ir .  The  ends  x  *  0  and  x  =  L  are  free  to 
deform. 

This  report  uses  the  following  four  solutions  (see  Appendix  B  and  Table  2) 
for  the  central  deflection  of  a  simply  supported  circular  cylindrical  shell  of 
radius  R,  length  L,  skin  thickness  h,  ana  material  with  Young's  Modulus  E  and 
poisson's  ratio  v  subject  to  a  point  load  P  (Figure  1).  The  deflection  (w)  is 
computed  under  the  load. 


ODQVIST  SOLUTION30 


V2-  >/2~[12(l-v2)]5/8  R3/4  L1/2 

2*  V57  h9/4 


n-1,3,5 
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where  Che  infinite  sum  can  be  calculated  as  1.68876336  correct  to  six  decimal 
places . 


BIJLAARD'S  SOLUTION36 


W  =  £  £  im2—'  v/-— —  L4  Z 
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for  m  =  0  and  n  =  1,3, 5, 7,...,  (4) 
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where 


y  *  mna/L 


A  =  1  +  Hr  *  B 


=  1  + 


fe(iT  .‘--fefeoo2 


F  - 


jl  +  (l-v)2l  /hf 

S  12(1-V)  |  ^Ry 


U  -  1  + 


and 


H  =  Aw8  +  4Bu6n2  +  6  Cy4n4  -  (8  -  2\>2)  y4n2  +  6Fy4  +  4u2  (n2-l)2n2 


+  n4  (n2  -  l)2 


YUAN-TING'S  SOLUTION25 
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A  = 


K1 

*1 

K2 

*2 

S 

*3 

K4 

*4 

nl 

n2 

n3 

\ 

n5 


(nx2  +  n22)(n12  +  n32) 


*l(n42  '  4Al\2)  +  442n4AlBl 


A1  (n i  -  n2n3)  -  B1n^^2  +  n3 


*2(n42  '  4Al2fil2)  ‘  4*ln4AlBl 


Bl^nl  ”  n2n3^  +  Aln 1 (n2  +  n3 


*3(„32  -  ^Gj2)  *  *V5C1G1 


Cl^nl  +  n2n3^  +  Glnl^n2  ”  n3 


\("32  -  4ClY>  *  4,3"5C1G1 


GlCn i  +  n2n3)  -  c1'i1^2  “  n3 


2  2  2  2 
A1  "  B1  "  C1  +  Gl 


2(A]B1  +  C1G1) 


2(A1B1  +  C1G1) 


2  2  2 
Ai  -B1  +n 


2  2  2 
C1  *  G!  +  n 


) 


) 


) 


) 


and 

,  C]i  =  real  parts  of  complex  roots 

Bl»  Gi  =  imaginary  parts  of  complex  roots  of  the  characteristic 
equation. 

Furthermore,  the  inextensional  solution-^  for  the  pinched  load  case 


p_ 

EL 
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This  solution  was  obtained  by  enforcing  the  inextensibility  condition  to  the 
three  strain  components  at  the  shell's  middle  surface.  The  assumed  tangential 
displacement  does  not  satisfy  the  boundary  conditions  at  the  free  ends.  We  note 
that  for  small  values  of  L,  Equation  (13)  is  of  the  same  form  as  the  "ring 
formula"  of  Reference  58  (p.  220).  However,  from  physical  grounds  the 
deflection  under  the  load  should  not  become  zero  as  L  increases,  but  have  a 
finite  value,  which  means  that  Equation  (13)  is  not  valid  for  large  L. 

To  complete  this  section,  we  note  that  the  case  of  an  infinitely  long  thin 
shell  subject  to  a  pinched  load  is  given  in  Buchwald  (p.  246  of  Reference  50) 
among  others.  The  deflections  under  the  load  and  at  ir / 2  away  from  it  are 
given  by: 


W(o) 


[3(l-v2)] 


3/4 


/  r\3/2  P_ 

Ih  I  Eh 
'  '  k=l 


(4k2-l) 


3/2 


TT 

wCJ) 


0.804168  (l-v2) 

4  .  3/4 

7  I3(l-v2)] 


(-l)kk 
,  9  ^/2 

(4k2-l) 


,  ox3/4 

*  -0.486281  ( 1-v2 ) 

Finally  a  FORTRAN  program  has  been  developed  for  the  VAX  11/780  computer  to 
obtain  the  maximum  central  deflection  of  a  simply  supported  cylindrical  shell 
subject  to  a  radial  point  load  by  the  methods  of  Odqvist,  3^  Bijlaard,3b 
Mizoguchi,3^  and  Yuan-Ting23  as  mentioned  previously.  In  this  program  the 
algorithm  for  any  cubic  equation,  based  on  the  standard  trigonometric 
method,33  was  developed  and  incorporated,  in  preference  to  approximate  or 
iterative  methods  such  as  those  shown  in  References  23,  45,  5,  66,  67,  and  78, 
which  in  the  precomputer  era  were  popular.  Furthermore,  the  inextensional 
solution  of  Timoshenko4  for  the  pinched  cylinder  case  at  midlength  with  free 
ends  has  also  been  included. 
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CHAPTER  4 

NUMERICAL  SOLUTIONS 


The  models,  which  were  studied  in  this  report,  were  discretized  by  using  a 
rectangular  grid  of  points,  when  the  elements  were  "plate  elements,"  or  a  curved 
rectangular  grid,  when  "shell  elements"  were  employed. 

POINT  LOAD  CASE 

For  STAGSC  the  discretization  was  carried  out  employing  "410"  plate 
quadrilateral  elements.  Their  characteristics  are  described  in  Reference  2. 

For  ABAQUS  the  modeling  was  carried  out  by  using  8-noded  "S8R"  shell  elements 
with  reduced  integration. ^  In  all  cases  only  half  of  the  shell  in  axial 
length  was  analyzed. 

There  were  two  models  used  with  both  STAGSC  and  ABAQUS  discretization.  The 
first  one  (90°)  was  employed  to  investigate  whether  a  simplification  in  halving 
the  peripheral  number  of  points  can  be  made  at  the  expense  of  accuracy.  The 
second  one  (180°)  is  theoretically  correct. 

The  discretization  for  STAGSC,  where  the  first  number  indicates  number  of 
axial  nodes  and  the  second  number  of  peripheral  nodes,  were: 

5  x  13  for  an  angle  of  90° 

5  x  25  for  an  angle  of  180° 

For  ABAQUS,  due  to  midside  nodes,  the  employed  grids  were: 

9  x  25  for  an  angle  of  90°. 

9  x  49  for  an  angle  of  180°. 

The  following  boundary  conditions  were  enforced: 

At  the  simply  supported  side  x=0  :  V  =  W  *  RU  »  0.* 

At  midlength  (x  *  L/2)  symmetry  :  U  *  RV  *  RW  *  0. 


*  U  =  Displacement  along  the  local  and  global  axial  axes  (coincident) 

V  »  Displacement  along  the  local  tangential  axis 

W  *  Displacement  along  the  local  radial  axis 

RU  *  Rotation  about  the  local  and  global  axial  axes  (coincident) 

RV  *  Rotation  about  the  local  tangential  axis 
RW  *  Rotation  about  the  local  radial  axis 
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At  <9  *  0°  (top  generatrix)  symmetry  :  V  =  RU  =  RW  =  0 . 
At  either  =  90°  or  180°  symmetry  :  V  =  RU  =  RW  =  0. 
All  rigid  body  modes  were  suppressed. 


PINCHED  LOAD  CASE 

For  STAGSC  the  modeling  was  carried  out  using  "410"  plate  elements,  while 
for  ABAQUS  through  "S8R"  shell  elements.  Again  only  half  of  the  length  of  the 
shell  axially  was  modeled.  Only  one  type  of  discretization  was  employed  in  this 
case.  It  is  given  below  for  both  STAGSC  and  ABAQUS. 

STAGSC  discretization:  5  x  13  for  an  angle  of  90°. 

ABAQUS  discretization:  9  x  25  for  an  angle  of  90°. 

The  following  boundary  conditions  were  applied: 

At  the  free  end  (x=0)  :  RV  =  RW  =  0 . 

At  raidlength  (x=L/2)  symmetry  :  U  =  RV  =  RW  =  0. 

At  «p  =  90°  (top  generatrix)  symmetry  :  V  *  RU  =  RW  ■  0. 

All  rigid  body  modes  were  suppressed. 
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CHAPTER  5 

COMPARISON  OF  NUMERICAL  AND  ANALYTICAL  SOLUTIONS 


The  objective  of  this  work  is,  while  trying  to  bring  together  some  of  the 
analytical  solutions  in  the  open  literature,  to  further  validate  the  two  finite 
element  programs  STAGS*  and  ABAQUS^  in  the  elastic  range.  For  this  purpose, 
a  series  of  ten  models  each  with  a  radius  R=100  inches,  thickness  h=1.0  inch, 
and  total  length  L  ranging  from  100.0  to  515.98  inches  has  been  analyzed.  Table 
3  displays  all  geometrical  material  properties.  Tables  4,  5,  8,  9,  10,  13,  ’5, 

^  and  17  pertain  to  the  point  load  case  (simple  supports,  V  =  W  =  RU  “  0).  Tables 

>  6,  7,  11,  12,  14,  and  16  contain  results  for  the  pinched  load  case  (continuous 

supports,  i.e.,  RV  =  RW  =  0). 


POINT  LOAD  CASE 

Table  4  summarizes  results  of  STAGS*  analyses  for  a  simply  supported 
shell  by  modeling  only  a  90°  segment  peripherally  and  half  length  axially. 
Although  we  know  that  this  is  not  exactly  rigorous,  at  least  for  static 
analysis,  we  would  like  to  examine  its  effect  on  the  maximum  central  deflection. 

Table  5  displays  STAGS*  results  by  modeling  the  correct  180°  segment 
(Figure  14).  It  can  be  seen  that,  as  L/R  increases,  the  two  answers  start 
deviating  by  as  much  as  4.5  percent. 

Table  9  displays  ABAQUS^  results  for  simply  supported  shells  subject  to  a 
point  force  radially  inwards  by  modeling  a  90°  segment  (Figure  15)  and  should  be 
compared  to  those  of  Table  4  (STAGS  values).  Table  10  contains  ABAQUS  results 
for  a  simply  supported  shell,  by  modeling  the  correct  180°  segment  and  should  be 
compared  to  Table  5  (STAGS  values).  The  STAGS  results  for  both  cases  are  lower 
(stiffer  model)  than  those  by  ABAQUS.  In  fact,  Figure  12  (Case  of  R/h  =  100) 
displays  this  feature  clearly.  We  may  observe  here  that  ABAQUS  employs 
Koiter-Sanders  theory  with  a  discretely  imposed  Kirchhoff  assumption  through  a 
penalty  method  formulation.  STAGS  on  the  other  hand  uses  a  first  order  shell 
theory  applicable  to  moderate  rotations.  The  differences  apparently  must  be 
attributed  to  the  differences  between  flat  (STAGS)  versus  curved  shell  elements 
(ABAQUS)  with  reduced  integration  and  not  so  much  to  the  formulation  since  we 
are  dealing  with  "linear  solutions." 

Tables  13  and  15  contain  analytical  results  for  the  range  of  radius/ 
thickness  (R/h)  of  100  to  207.36  and  radius/ length  (R/L)  from  0.0969  to  1.0. 
Table  17  covers  the  range  of  R/h  from  200  to  414.7  and  R/L  from  0.1938  to  1.0. 
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Figures  3  to  6,  and  8  are  plots  of  results  of  Tables  13  and  17, 
respectively,  in  terms  of  EWH/P  versus  R/L  for  various  R/h  ratios.  Figures  9 
and  10  are  extensions  of  Figures  5  and  6,  respectively,  over  R/h  ratios  from  200 
to  414.72.  Finally  Figure  12,  which  applies  to  the  case  of  R/h=100,  summarizes 
the  various  analytical  and  numerical  answers  for  the  simply  supported  case. 
Generally,  computations  by  STAGS  yielded  a  stiffer  model.  The  apparent 
increasing  disagreement  for  ABAQUS  results  (Figure  12)  from  the  analytical 
solutions  as  L/R  increases  are  due  to  discretization.  In  all  cases,  the  same 
mesh  was  used;  it  must  be  stressed,  however,  that  the  ABAQUS  elements  with  their 
midside  nodes  increased  the  problem  size  considerably,  while  increasing  accuracy. 

PINCHED  LOAD  CASE 

Table  6  displays  pinched  load  results  by  STAGS.  Tables  7  and  8  are  two 
convergence  studies  for  both  the  pinched  cylinder  and  the  simply  supported 
cylinder  cases.  The  two  answers  do  not  vary  significantly. 

Table  11  reports  ABAQUS  results  for  the  pinched  cylinder  problem  and  must 
be  compared  to  Table  6  (STAGS). 

Tables  12  and  14  give  the  analytical  results  for  R/h  from  100.0  to  207.36, 
and  R/L  from  0.0969  to  1.000.  Table  16  extends  the  range  for  R/h  from  200.0  to 
414.7  and  R/L  from  0.1938  to  1.000.  Figures  7  and  11  show  curves  corresponding 
to  the  pinched  cylinder  problem,  as  given  by  Tables  12  and  16  respectively. 
Finally  Figure  13,  which  applies  to  the  case  of  R/h=100,  summarizes  the  various 
analytical  and  numerical  answers  for  the  pinched  case.  For  this  case  (Figure 
13)  both  STAGS  and  ABAQUS  produced  larger  displacements  than  the  analytical 
solution.  Here  STAGS  was  closer  to  theory  than  ABAQUS. 

Based  on  comparisons  with  the  computed  analytical  solutions,  ABAQUS 
performed  better  for  the  simply  supported  case,  while  STAGS  performed  better  for 
the  pinched  cylinder  problem.  Naturally  it  would  have  been  much  better  if  the 
previous  computational  and  analytical  answers  were  verified  by  experiments, 
which  is  the  ultimate  check.  What  has  emerged  from  the  above  analyses  is  the 
totally  different  behavior  for  the  two  types  of  support  conditions  and  loadings 
(simply  supported  versus  pinched  cylinders)  attributed  to  the  membrane  action 
even  in  the  elastic  range.  The  implication  is  this,  although  the  behavior  is 
different,  for  the  same  R/h  as  L/R  increases,  it  "appears"  that  both  solutions 
tend  to  the  same  value.  However,  the  point  load  case  approaches  from  below  some 
asymptotic  value,  while  the  pinched  load  case,  as  can  be  seen  from  Equation  (13), 
approaches  zero  asymptotically  from  above,  which  is  not  physically  correct. 
Seide's  solution^  addresses  the  infinitely  loaded  cylinder  subject  to  pinched 
radial  loads.  From  the  practical  point  of  view,  however,  simplification  of  a 
problem  to  continuous  supports  from  simple  supports  cannot  occur,  except  in  the 
case  of  a  pipeline,  which  is  indeed  a  long  shell. 
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CONCLUSION 

This  study  brings  together  existing  analytical  solutions  for  the  cases  of 
(1)  a  simply  supported  linear  elastic  circular  cylindrical  shell  subject  to  a 
radially  inwards  concentrated  force  at  midspan  and  (2)  a  linear  elastic  pinched 
cylindrical  shell  at  midspan  with  free  ends.  It  compares  the  analytical 
solutions  with  finite  element  computations  by  both  STAGS ^  and  ABAQUS^ 
general  purpose  computer  programs.  Agreement  is  fairly  good  for  both  codes. 
While  the  ABAQUS  answers  are  nearer  to  the  analytical  answers,  it  must  be  added 
that  computations  by  STAGS  were  substantially  faster.  Furthermore,  for  static 
analyses,  where  the  characteristic  length  circumferentially  (such  as  2R  for  a 
180°  modeling,  versus  R  for  a  90°  segment)  does  not  affect  the  calculation  (for 
frequency  or  shock  analyses  this  is  not  true) ,  it  has  been  shown  by  these  simple 
calculations,  that  the  effect  of  the  load  is  indeed  very  localized  (as 
demonstrated  by  the  effect  of  the  characteristic  roots  of  the  homogeneous 
equation). 5  As  far  as  maximum  displacement  (in  the  elastic  range)  at  the  point 
of  application  is  concerned,  it  is  sufficient  to  model  a  90°  segment  only.  This 
is  of  importance  to  structural  designers,  as  it  reduces  computational  effort 
substantially. 
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FIGURE  5.  NONDIMENSION AL  PARAMETER  EHW/P  FOR  A  SIMPLY  SUPPORTED  CYLINDRICAL  SHELL  SUBJECT 
TO  A  RADIALLY  INWARDS  CENTRALLY  LOCATED  CONCENTRATED  FORCE  P  (BY  MI20GUCHI 
SOLUTION)  PLOTTED  VERSUS  R/L  FOR  VARIOUS  R/H  RATIOS 


...  .IS  .2,  3«  •«>  ^  98  •9#  '•'8  '-lt 


FIGURE  6.  NONDIMENSION AL  PARAMETER  EHW/P  FOR  A  SIMPLY  SUPPORTED  CYLINDRICAL  SHELL  SUBJECT 
TO  A  RADIALLY  INWARDS  CENTRALLY  LOCATED  CONCENTRATED  FORCE  P  (BY  YUAN-TING 
SOLUTION)  PLOTTED  VERSUS  R/L  FOR  VARIOUS  R/H  RATIOS 


» 


nriafus/iENbTH 


FIGURE  10.  NONDIMENSIONAL  PARAMETER  EHW/P  FOR  A  SIMPLY  SUPPORTED  CYLINDRICAL  SHELL 
SUBJECT  TO  A  RADIALLY  INWARDS  CENTRALLY  LOCATED  CONCENTRATED  FORCE  P  (BY 
YUAN-TING  SOLUTION)  PLOTTED  VERSUS  R/L  FOR  VARIOUS  R/H  RATIOS 
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H 


FIGURE  12.  SUPERIMPOSED  PLOTS  OF  ANALYTICAL  SOLUTIONS  BY  FOUR  METHODS  AND  FOUR  FINITE 
ELEMENT  DISCRETIZATIONS  OF  EHW/P  VERSUS  L/R  FOR  R/H  =  100 


FOR  VARIOUS  METHODS 
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FIGURE  14. 


HOCEL  SCRLE  =  3.M5!0E-»H.  ORIENT.  -  31.  1J. 

5CLJTICK  ECRLC  •  3.77-OCC-tm 


SIMPLY  SUPPORTED  CYLINDER  (  9  -  180°)  SUBJECTED  TO  RADIAL  CONCENTRATED  LOAD 
(BY  STAGS) 
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TABLE  1.  GOVERNING  DIFFERENTIAL  EQUATIONS  USED  FOR  ANALYTICAL  SOLUTIONS* 
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TIMOSHENKO^*  Energy  Argument  Through  Principle  of  Virtual  Displacements 


*  For  definition  of  employed  variables  see  Nomenclature 


AUTHOR  REF  EQUATION  TYPE  PROBLEM  TYPE  BOUNDARY  CONDITIONS  USED  METHOD 
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TABLE  3.  GEOMETRICAL  AND  MATERIAL  PROPERTIES  OF  ANALYZED  MODELS 


I 


M  - 
h 
fc 
% 

% 

*  . 

MODEL 

RADIUS 

R 

THICKNESS 

h 

ANALYZED 

LENGTH 

L 

TOTAL 

LENGTH 

L'j* 

YOUNG’S 

MODULUS 

E 

POISSONS  * 

RATIO  > 

NO. 

(in. ) 

(in. ) 

(in.  ) 

(in. )  R/h 

Lt/R 

(Lb/ in-) 

V  1 

50.000 

100.00 

100 

1.0000 

30 

X 

10u 

0.300 

60.000 

120.00 

100 

1.2000 

30 

X 

106 

0.300 

72.000 

144.00 

100 

1.4400 

30 

X 

io6 

0.300 

86.400 

172.80 

100 

1.7280 

30 

X 

106 

0.300 

103.680 

207.36 

100 

2.0736 

30 

X 

io6 

0.300 

124. 415 

248.83 

100 

2.4883 

30 

X 

io6 

0.300 

149.300 

298.60 

100 

2.9860 

30 

X 

iob 

0.300 

179.160 

358.32 

100 

3.5832 

30 

X 

106 

0.300 

214.990 

429.98 

100 

4.2998 

30 

X 

io6 

0.300 

257.990 

515.98 

100 

5.1598 

30 

X 

106 

0.300 
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TABLE  4.  RESULTS  BY  STAGSC  COMPUTER  PROGRAM  FOR  A  SIMPLY  SUPPORTED 
SHELL  (V=W=RU=0)  SUBJECT  TO  A  CENTRALLY 
LOCATED  INWARD  CONCENTRATED  FORCE 


14.489 

-0.0002 

103.526 

13.176 

-0.0002 

113.843 

11.988 

-0.0002 

125.125 

10.927 

-0.0002 

137.274 

9.9609 

-0.0002 

150.588 

9.0331 

-0.0002 

166.055 

8.1473 

-0.0002 

184.110 

7.3972 

-0.0002 

202.779 

6.8559 

-0.0002 

218.789 

6.4700 

-0.0002 

231.839 
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TABLE  5.  RESULTS  BY  STAGSC  COMPUTER  PROGRAM  FOR  A  SIMPLY  SUPPORTED 
SHELL  (V=W=RU=0)  SUBJECT  TO  A  CENTRALLY 
LOCATED  INWARD  CONCENTRATED  FORCE 


NO.  OF  NODES  NO.  OF  NODES 


MODEL 

NO. 


AXIALLY 
INCLUDING 
MIDSIDE 
NODES 
(IF  ANY) 


PERIPHERALLY 
INCLUDING 
MIDSIDE 
NODES 
(IF  ANY) 


ANALYZED 
SUBTENDED 
ANGLE 
(  °) 


MIDSIDE 

NODES 


REACTION  RADIAL 

ON  QUARTER  INWARD 

MODEL  DISPLACE 

P/4  MENT 

(LB)  (IN) 


1 

5 

25 

180° 

NO 

-14.489* 

-0.0002 

103.526 

2 

5 

25 

180° 

NO 

-13.180 

-0.0002 

113.808 

3 

5 

25 

180° 

NO 

-11.992 

-0.0002 

125.083 

4 

5 

25 

180° 

NO 

-10.911 

-0.0002 

137.475 

3 

3 

25 

180° 

NO 

-  9.930 

-0.0002 

151.057 

6 

5 

25 

180° 

NO 

-  9.0428 

-0.0002 

165.877 

7 

5 

25 

180° 

NO 

-  8.2397 

-0.0002 

182.045 

8 

5 

25 

1— * 

00 

o 

o 

NO 

-  7.5112 

-0.0002 

199.701 

9 

5 

25 

180° 

NO 

-  6.8466 

-0.0002 

219.086 

10 

5 

25 

180° 

NO 

-  6.2267 

-0.0002 

240.898 

*A11  values 
in  element 
-13.238. 

were 

410) 

obtained  with 
in  N-l  cards. 

IPENL  =  0 
If  IPENL 

(penalty  function  on  fourth  order 
=  1  the  first  entry  for  P/4  become 

terms 

s 
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TABLE  6.  RESULTS  BY  STAGSC  COMPUTER  PROGRAM  FOR  A  PINCHED 
SHELL  (BOUNDARY  CONDITIONS  AT  ENDS  RV=RW=0) 


MODEL 

NO. 

NO.  OF  NODES 
AXIALLY 
INCLUDING 
MIDSIDE 

NODES 
(IF  ANY) 

NO.  OF  NODES 
PERIPHERALLY 
INCLUDING 
MIDSIDE 

NODES 
(IF  ANY) 

ANALYZED 

SUBTENDED 

ANGLE 

(°) 

MIDSIDE 

NODES 

REACTION 

ON  QUARTER 
MODEL 

P/4 

(LB) 

W 

RADIAL 

INWARD 

DISPLACE¬ 

MENT 

(IN) 

Ehw 

P 

1 

5 

13 

90° 

NO 

-0.18425* 

-0.0002 

8141.112 

2 

5 

13 

90° 

NO 

-0.22068 

-0.0002 

(797.172 

3 

5 

13 

90° 

NO 

-0.26414 

-0.0002 

5678.806 

4 

5 

13 

90° 

NO 

-0.31592 

-0.0002 

4748.037 

5 

5 

13 

90° 

NO 

-0.37746 

-0.0002 

3973.931 

6 

5 

13 

90° 

NO 

-0.45039 

-0.0002 

.'330.446 

7 

5 

13 

90° 

NO 

-0.53655 

-0.0002 

2795.638 

8 

5 

13 

90° 

NO 

-0.63789 

-0.0002 

2351.502 

9 

5 

13 

90° 

NO 

-0.75636 

-0.0002 

3983.182 

10 

5 

13 

90° 

NO 

-0.89353 

-0.0u02 

3678.734 

*AiI  values  were  obtained  with  IPENL  =  0  (penalty  function  on  fourth  order  terms 
in  element  410)  in  N-l  cards.  If  IPENL  =  1  the  first  entry  for  P/4  becomes 
-0.18404. 
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TABLE  7.  CONVERGENCE  STUDY  BY  STAGSC  COMPUTER  PROGRAM 
OF  MODEL  NO.  1  WITH  CONTINUITY  END  CONDITIONS 


MODEL 

NO. 

NO.  OF  NODES 
AXIALLY 
INCLUDING 
MIDSIDE 
NODES 
(IF  ANY) 

NO.  OF  NODES 
PERIPHERALLY 
INCLUDING 
MIDSIDE 
NODES 
(IF  ANY) 

ANALYZED 

SUBTENDED 

ANGLE 

(°) 

MIDSIDE 

NODES 

REACTION 

ON  QUARTER 
MODEL 

P/4 

(LB) 

W 

RADIAL 

INWARD 

DISPLACE 

MENT 

(IN) 

1 

5 

13 

90° 

NO 

-0.18425 

-0.0002 

1 

9 

25 

90° 

NO 

-0.18361 

-0.0002 

l 

19 

51 

90° 

NO 

-0.18343 

-0.0002 

Ehw 

P 

8141.112 

8169.489 

8177.506 
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TABLE  8.  CONVERGENCE  STUDY  BY  STAGSC  COMPUTER  PROGRAM 
OF  MODEL  NO.  1  WITH  SIMPLY  SUPPORTED  ENDS 


NO.  OF  NODES 

NO.  OF  NODES 

W 

AXIALLY 

PERIPHERALLY 

REACTION 

RADIAL 

INCLUDING 

INCLUDING 

ANALYZED 

ON  QUARTER 

INWARD 

MIDSj.DE 

MIDSIDE 

SUBTENDED 

MODEL 

DISPLACE¬ 

MODEL 

NODES 

NODES 

ANGLE 

MIDSIDE 

P/4 

MENT 

NO. 

(IF  ANY) 

(IF  ANY) 

(°) 

NODES 

(LB) 

(IN) 

^3  3" 


1  jf  ww  iw  w  11  w  w  w  w  juiwj 

1 
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TABLE  9.  RESULTS  BY  ABAQUS  COMPUTER  PROGRAM  FOR  A  SIMPLY 

SUPPORTED  SHELL  (V=W=RU=0)  SUBJECT  TO  A  CENTRALLY 
LOCATED  INWARD  CONCENTRATED  FORCE 


NO.  OF  NODES 

NO.  OF  NODES 

W 

AXIALLY 

PERIPHERALLY 

REACTION 

RADIAL 

INCLUDING 

INCLUDING 

ANALYZED 

ON  QUARTER 

INWARD 

MIDSIDE 

MIDSIDE 

SUBTENDED 

MODEL 

DISPLACE- 

MODEL 

NODES 

NODES 

ANGLE 

MIDSIDE 

P/4 

MENT 

NO. 

(IF  ANY) 

(IF  ANY) 

(°) 

NODES 

(LB) 

(IN) 
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TABLE  10.  RESULTS  BY  ABAQUS  COMPUTER  PROGRAM  FOR  A  SIMPLY 

SUPPORTED  SHELL  (V=W=RU=0)  SUBJECT  TO  A  CENTRALLY  LOCATED 
INWARD  CONCENTRATED  FORCE 


NO.  OF  NODES  NO.  OF  NODES 


MODEL 

NO. 


AXIALLY 
INCLUDING 
MIDSIDE 
NODES 
(IF  ANY) 


PERIPHERALLY 
INCLUDING 
MIDSIDE 
NODES 
(IF  ANY) 


ANALYZED 

SUBTENDED 

ANGLE 

(°) 


MIDSIDE 

NODES 


REACTION 
ON  QUARTER 
MODEL 
P/4 
(LB) 


W 

RADIAL 

INWARD 

DISPLACE¬ 

MENT 

(IN) 


Ehw 

P 


i'|  VL/iA  «**  .  a  i  »-«j 


.1 .  .1  Uq 
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TABLE  11.  RESULTS  BY  ABAQUS  COMPUTER  PROGRAM  FOR  A  PINCHED 
SHELL  (BOUNDARY  CONDITIONS  AT  ENDS  RV=RW=0) 


NO.  OF  NODES  NO.  OF  NODES 
AXIALLY  PERIPHERALLY 


MODEL 

NO. 


INCLUDING 
MIDSIDE 
NODES 
(IF  ANY) 


INCLUDING 
MIDSIDE 
NODES 
(IF  ANY) 


ANALYZED 

SUBTENDED 

ANGLE 

(°) 


MIDSIDE 

NODES 


W 

REACTION  RADIAL 

ON  QUARTER  INWARD 

MODEL  DISPLACE- 
P/4  MENT 

(LB)  (IN) 
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TABLE  12. 


PARAMETER  EWH/P  FOR  VARIOUS  R/H  AND  R/L  RATIOS  FOR  A  PINCHED  CYLINDRICAL  SHELL  BY 
TIMOSHENKO'S  INEXTENSIONAL  SOLUTION 


EAHAH/P 

r/h  =  a 

RADIUS 

THICKNESS 

LENGTH 

R/L  INEXTENSIONAL 

IBB. 8BB 

1.  BBB 

IBB  BB 

B. 1B0B00E+01 

B. 812330E+B4 

100000 

1.  BBB 

1SBBB 

B  833333E+B0 

B  474942E+B4 

100  000 

1  BBB 

144  BB 

B. 494444E+BB 

8  544118E+B4 

IBB. BBB 

1  BBB 

172  88 

B.  5787B4E4BB 

B. 47BB98E+B4 

lea  aae 

1  BBB 

2B7. 34 

8. *8225 3E+BB 

B  391749E+B4 

IBB. BBB 

1  BBB 

248.  83 

B.  401878E+0B 

B. 324457E+B4 

IBB  BBB 

1  BBB 

298  48 

8.  334898E+BB 

8. 272B48E+B4 

IBB. BBB 

1  BBB 

358 . 32 

6. 279082E+00 

B. 226706E+B4 

IBB  BBB 

1  BOB 

429.98 

B. 232548E+BB 

0  188922E4B4 

IBB. BBB 

1 .  BBB 

S 15 . 98 

B.  19  3887E4BB 

B. 157435E+B4 

R/H  =  a 

120 

00® 

1 

.  00® 

120 

00 

8. 100800E+01 

0 

1 16976E+05 

120 

000 

1 

000 

144 

00 

8  833333E+00 

0 

974796E+04 

120 

000 

1 

00® 

172 

a® 

B  494444E400 

0. 

8 12  330E  +  04 

120 

000 

1 

000 

207 

36 

B  575784E4BB 

0 

*769422*04 

120 

000 

1 

000 

248 

83 

8  48225  3E  400 

0 

5641 18E+04 

120. 

000 

1 

000 

2*6 

*0 

8  *81S79£f88 

0 

4700922*04 

120 

000 

1 

00® 

358 

32 

B  334898E+BB 

0 

391749E+04 

120 

00® 

1 

000 

429 

98 

B  279B82E4BB 

0 

32645  7E  +  04 

120 

000 

1 

000 

515 

98 

8  232568E480 

0 

272048E  *04 

120 

000 

1  ■ 

000 

619 

17 

3. i933e7E+30 

0 

226706E+04 

*/b 

1  =  0 

1440SCE+03 

3  44 

000 

1 

000 

i** 

00 

a  wBsaaaE+ai 

0 

1 6S445E+ 05 

144 

000 

1 

000 

172 

80 

8. 333333E+88 

0 

140371E+0S 

144 

000 

1 

000 

287 

26 

8. 694444E+00 

0 

114974E+05 

144 

000 

1 

000 

2*8 

83 

8.  !i7B78*E+B0 

0 

974794E+04 

144 

000 

1 

000 

298 

*0 

a  *8225  3E+00 

0 

812330E+04 

144 

000 

1 

000 

358 

32 

8  401878E+00 

0 

474942E+04 

1  44 

000 

1. 

000 

*29 

98 

8  33*898E+00 

0 

5  6*1 1 8E  +  0* 

144 

000 

1 

000 

SIS 

98 

0 . 279082E+00 

0 

470098E+04 

144 

000 

1 

000 

819 

17 

0  23254SE+00 

0 

391749E+04 

144 

000 

1 

000 

7*3 

01 

0. 193807E+88 

0 

324457E+04 

R/ 

H  =  0 

172388E+03 

172 

800 

1 

000 

172 

60 

0  l 00000E+0 1 

0 

242560E+0S 

172 

600 

1 

000 

207 

36 

0  333333E+00 

0 

202134E+05 

172 

300 

1 

000 

248 

33 

0 . 694444E+A0 

0. 

16S445E+05 

172. 

800 

1 

000 

293 

60 

0  378704E+08 

0 

140371E+05 

172 

800 

1 

00® 

3S8 . 

32 

0  *8225  3E+00 

0 

1 14974E+05 

172 

300 

1 

000 

429 

98 

0  10187SE+00 

0 

974794E+04 

172 

300 

1 

000 

515 

98 

0  334898E+00 

0 

812330E+0* 

172 

800 

1 

000 

619 

17 

0.  279082E+00 

0 

474942E+04 

172 

30® 

1 

000 

743 

01 

0  232548E+00 

0 

54*1 18E+04 

172. 

80® 

1 

000 

391 

61 

0.  193807E+00 

0 

470098E+0* 

B/H  =  0 

207340E+03 

207  340 

1 

080 

207. 

36 

0. 100000E+01 

0 

349287E+05 

207  340 

1 

000 

2*3 

83 

0  333333E+00 

0 

291073E+05 

207  340 

1 

000 

298 

60 

0. 494444E+00 

0 

2425  40E+05 

20 7  340 

1 

000 

358 

32 

0  5  78704E+00 

0 

202134E+0S 

207  340 

1 

000 

*29 

98 

8  48225 3E+00 

0 

1 48445E+05 

207  340 

1 

000 

515 

98 

0. 401878E+00 

0 

140371E+0S 

£07. 340 

1 

000 

419 

17 

0  334898E+00 

0 

1 1 4974E+05 

207. 340 

1 

000 

7*3 

01 

0  279082E+00 

0 

974794E  f04 

207  340 

1 

000 

891 

61 

0  232S48E+00 

0 

8123 j0E+04 

207  340 

1 

000 

1049 

93 

0  193S07E+00 

0 

474942E+04 

42 
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TABLE  13.  PARAMETER  EWH/P  FOR  VARIOUS  R/H  AND  R/L  RATIOS  FOR  A  SIMPLY  SUPPORTED  CYLINDRI¬ 
CAL  SHELL  SUBJECT  TO  A  CENTRALLY  LOCATED  INWARD  POINT  LOAD  BY  FOUR  ANALYTICAL 
METHODS 

£*W*H/P 


R'rt  s  0  1O00O0E+03 

RADIUS  THICKNESS  LENGTH  R/L  0D9VIST  BI7LAARO  HIZOOUCHI  YUAN— TING 


100  000  1.000  100.00  0  100000E+01  0. 115627E+03  0.11 IS95E+03  0 . 1 10917E+03  0  U2321E+03 
100  000  1.000  120.00  0  S33333E+00  0.  126663E+03  0  1231S6E403  0.1224S9E+03  0  123925E+S3 
100.000  1.000  1 44 . 00  0.S94444E+00  0  13B7S2E+03  0  135860E+03  0  13S141E+03  8  136630E+03 
100.000  1 .000  172.80  0.  578704E+00  0.151995E  +  03  0  149839E+03  0.149099E+03  0  . 150369E+03 
100.000  1.000  207.34  0.4822S3E+00  0.166502E+03  0  16S243E+03  0  164479E+03  0  16S543E+03 
100008  1.000  248.83  0. 401873E+00  0. 132394E+03  0. 192243E+03  O.1814S0E703  0.133277E+03 
100  000  1.000  293.60  0. 33489SE+00  0  199303E+03  0 . 20105 9E403  0  200222E+03  0.204024E+03 
100  000  1  800  358. 32  0  279082E+00  0  218873E+03  a  221379E+03  0  220976E+03  0  22S798E+03 
100  000  1.000  429  93  0  23254SE+00  0 . 239763E+03  0  244962E+03  0  243952E+03  0  24S4O9E+03 
100  000  1.000  SIS. 98  0. 193807E+00  0  262643E+03  0  27097SE+03  0  269789E403  0 . 26279 3E+03 


R/H  =  0  120O00E+03 


120.000  1.000  120  00  0.  100000E+01  0.  14S223E  +  03  0.J397S7E+03 
128  800  1.000  144  00  0. 833333E+00  0  1S9083E+03  0  IS4217E+03 
120  000  1.800  172.30  0  O94444E+00  0.174267E+03  0  170101E+03 
120  000  1  000  207  36  0. S78704E+00  0  190900E+O3  0. 137S67E+03 
120  000  1  000  248  83  0. 48225 3ET00  0  209121E+03  0 . 20t8O2E+03 
120  000  1.000  298  60  0. 401878E+00  0  2290B0E+03  0 . 228007E+03 
120  800  1.000  358.32  0. 334S98E+00  0  2S0945E+03  0  251431E+03 
120  0O3  1  000  429  93  0 . 279082E+00  0. 274896E+03  0  277367E+03 
120  000  1  000  515  98  0. 232S62E+0O  0.301134E+03  0.38e011E+03 
120.000  1  000  619.17  0. 193807E+00  0  329875E+03  0.337954E403 


R/H  =  0. 144000E+03 


144  000  1.000 
144  000  1 . 000 
144  000  1.000 
144  000  1.000 
144  000  1 . 000 
144  800  1.000 
144  000  1.000 
144.000  1.000 
144  300  1. 000 
144  000  1.000 


R/H  =  0.172300E+03 

172  800  1. 000 
172  800  1.000 
1  72  800  1 . 000 
172  300  1.000 
172  300  1.000 
172  300  1.000 
172  S00  1 . 800 
172  800  1  000 
172  900  1  000 
172  300  1. 000 


R/H  =  O. 207360E+03 


207  360 

1. 000 

207  36 

0. 100000E+01 

0 

287716E+03 

0. 273444E+03 

0 

E72779E+03 

0 

E8E393E+03 

207  360 

1 . 000 

243  83 

0. 333333E+00 

0 

315 177E+03 

0.  30i",'?4E‘p03 

0 

301106E403 

0 

3107S4E403 

207  360 

1  000 

298. 60 

0. 694444E+00 

0 

345259E+03 

0  33c900E+03 

0. 

332191E+03 

0 

341*35  *£+03 

207  360 

1. 000 

353.  32 

0. 573704E400 

0 

378213E+03 

0  367037E+03 

0 

366327E*03 

0 

376235Er03 

207  360 

1 . 000 

429. 98 

0. 48225 3E+00 

0 

4143I1E+03 

0. 40459»E+03 

0 

4C3237E 403 

0 

4141»£E403 

207. 360 

1 . 000 

SIS  98 

0.  401 878E+00 

0 

45  33S5E+03 

0.  445  S70E+03 

0 

445OS1E+03 

0 

454330E403 

207  360 

1  000 

619  17 

0  334898E+00 

0 

497173E+03 

0  491279E+03 

0 

49O45»E+03 

0 

4*«204Et-e? 

207  360 

1  000 

743  01 

0. 2790S2E+00 

0 

544626E403 

0.  S 4 1  30O>E+0 3 

0 

S40392E+03 

0. 

551448E403 

207  360 

1. 000 

891  61 

0.  232563E4O0 

0 

S96608E+03 

0  5«646*E403 

0 

5  95  4  30E  +  0  3 

0 

4l3^*tE403 

207. 360 

1. 000 

1069  93 

0  1 9  3807E+00 

0 

65  355 1E+0  3 

0  6S  7 2 16E+0  3 

9 

6SS993E+03 

a 

68045 7E403 

172  80  0 . 100000E+01  0.229080E+03  0  2187B7E+03  0  218071E+03  0.224362E+03 
207  36  0. 333333E+00  0  25094SE+03  0  241433E+03  0  24O692E+03  0  24699SE+03 
24883  0.694444E+00  0274896E+03  0266288E+03  026S526E+03  0. 271900E+03 

298.60  0.S78704E+00  0.301134E+03  0.293593E+03  0. 292S10E+03  0  29944SE+03 
358  .32  0. 48225  3E400  0. 329875E+03  0  323618E+03  0  322ei0E+03  0  329406E+03 
429  98  0. 401878E+00  0.361360E+03  0  J56667E+03  0  355827E+03  0  361597E+03 
SIS  98  0  334898E+00  0 . 39585 1E403  0.  393072E+03  0  392185E+03  0. 398034E+03 
619.17  0  279082E+00  0  433633E+03  e.43222°E+03  0.432269E+03  0  441833E+03 
743  01  0  2325 6SE+00  0.475021E+03  0 .  4776 02E+03  0.476526E+03  0  491784E+03 

391.61  0. 193807E+00  0.S203S9E+03  0  S26457E+03  0  525198E+03  0  540260E+03 


14400  0  10000OE+01  0 . 182394E+03  0.174922E+03  0  174196E+03  0. 17320OE+03 
172  80  0-  B33333E+00  0.199803E+03  0.  193015E+03  0.  19226SE+03  0.  196292E1-03 
207  36  0.694444E+00  0.218873E+03  0.212881E+03  0.212109E+03  0  216262E+03 
248.83  0.S7B704E+O0  0.239763E+03  0. 234716E+03  0  233923E403  0 . 236195Et*3 
298.60  0. 48225 3E+00  0. 262648E+03  8.253742E+03  0257925E+03  0261316E+03 
358.32  0  401873E+00  0.287716E+03  0. 285211E+03  0  284363E+03  0.287726E+03 
429.98  0. 334898E400  0.315177E+03  0. 3144O1E603  0. 313S07E+03  0. 318148E+03 
SIS  98  0- 279082E+00  0  3452S9E403  8.346674E+03  0  34S709E+03  0  354205E403 
619  17  0  232S6SE+00  0. 378213E+03  S. 332323E+B3  C.3312S1E+33  0.392322E+03 
743.01  0.193807E+00  0.414311E+03  0.421721E403  0  420465E+03  0 . 426385E403 


0  139047E+03 
0  IS  3435E+03 
0  169347E+03 
0  186793E+03 
0. 206003E+03 
0. 227179E+03 
0. 250561E+O3 
0  27*425  E  +  0  3 
0  30496OE403 
0. 336727E+03 


0  141 494E+03 
0  155985E+03 
0  171955E+03 
0  139312E+03 
0  208082E+03 
0. 22S38SE+03 
0. 2S4767E+03 
0  283416E+03 
0  311190E+03 
0. 334995E+03 
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TABLE  14.  PARAMETER  EWH/P  FOR  VARIOUS  R/H  AND  R/L  RATIOS  FOR  A  PINCHED  CYLINDRICAL 
SHELL  BY  TIMOSHENKO'S  INEXTENSIONAL  SOLUTION 

E9U9H/P 


r/h  s  a. 100000E+03 


RADIUS 

thickness 

LENGTH 

R/L 

INEXTENSIONAL 

100  000 

i  000 

200. 00 

0.  S00000E+00 

0  406165E+0* 

100. 000 

l  000 

2*0  00 

0. 414447E+00 

0  338471E+0* 

100  000 

i.  000 

288. 00 

0. 347222E+00 

0  28205 9E+0* 

100. 000 

l.  000 

3*5 . AO 

0. 2893S2E+00 

0. 235049E+0* 

100. 000 

l  000 

414.72 

0. 2411278*00 

0. 19S874E+0* 

100. 000 

l.  000 

497.44 

0. 200939E+00 

0  163229E+0* 

100  000 

1  000 

597.  20 

0.  147449E+00 

0  1 36024E+0* 

100. 000 

1  000 

714. 44 

0.  1 395  41E+00 

0.  11335  3E  +  04 

100. 000 

1  000 

859. 94 

0.  U6234E+00 

0. 944410E+03 

100. 000 

1  000 

1031. 94 

0. 949033E-01 

0.  787175E+03 

r/h  s  a 

120000E+O3 

120  000 

1  000 

2*0  00 

0  500000E+00 

0.  S84878E+0* 

120. 000 

1  000 

288  00 

0.  414447E+00 

0  487398E+0* 

120. 000 

1  000 

3*5  40 

0.  347222E+00 

0  406165E+0* 

120. 000 

1  000 

*1*  72 

0  289352E+00 

0. 338471E+0* 

120. 000 

1  000 

*97. 44 

0. 241 127E+O0 

0. 28205 9E+0* 

120. 000 

1  000 

S97. 20 

0  200939E+00 

0  235049E+04 

120  000 

1  000 

714. 44 

0  147449E+00 

0  195874E+0* 

120. 000 

1  000 

859 . 94 

0.  1 39541E+00 

0  143229E+0* 

120. 000 

1. 000 

1031. 94 

0.  114284ETO0 

0. 13402AE+0* 

120  000 

1.000 

1238. 35 

0  949033E-01 

0.  1133S3E+0* 

R/H  s  0. 

14A000E+03 

144.000 

1  000 

288  00 

0  500O00E+0O 

0  8*222*E+0* 

144  000 

1  000 

345  40 

•  •6667E+00 

0  701BS3E+0* 

144  000 

1. 000 

414  72 

0  3*7222E+0O 

0  584878E+0* 

144.000 

1  000 

497  64 

0  289352E+00 

0  487398E+04 

144  000 

1  000 

597  20 

0.  241 127E+00 

0  40616SE+0* 

144  000 

1.000 

716  64 

0. 200939E+O0 

0  338471E+0* 

1*4.000 

1  000 

859  96 

0. 167449ET00 

0  28205 9E+0* 

144. 000 

1. 000 

1031  96 

0.  139541E+00 

0  23S049E+04 

144  000 

1  000 

1238  35 

0.  U62S4E+00 

0  1 95  874E+0* 

1*4  000 

1.000 

1486  02 

0 . 969033E-01 

0  163229E+0* 

R/H  3  0 

172800E+03 

172. 800 

1  800 

345  60 

0. S0O00OE+00 

9 

1212B0ET0S 

172  800 

1. 000 

414  72 

0  416667E+00 

9 

10106 7E+05 

172.800 

1  000 

497  66 

0. 347222E+00 

9 

84£22*E*04 

172.800 

1  000 

597  20 

0  289  35  3E+00 

0. 

70185 3E+0* 

172  300 

1  000 

716. 64 

0  241127E+00 

9. 

584878E+04 

172. 800 

1  000 

859  96 

0. 200939E+00 

9 

487398E+0* 

172.800 

1  000 

1031  96 

0. 167449E+00 

9 

406165E+0* 

172.800 

1  000 

1238  35 

0. 139541E+00 

9. 

338471E+0* 

172  800 

1  000 

1486  02 

0  116284E+00 

0 

2820S9E+0* 

172. 800 

1.000 

1783  22 

0  969033E-01 

9 

23S049E+04 

R/H  s  0. 

207360E+03 

207  360 

1  000 

*1*  72 

0  SO8000E+00 

0. 17*6**E+05 

207  340 

1  000 

*97  66 

0. 416667E+00 

0  14S5  36E605 

207  340 

1 . 000 

597  20 

0.  347222E40O 

0  121280E+05 

207  360 

1  000 

716  6* 

0.  289352E+00 

0. 101067E+05 

207  340 

1  000 

8S9.  96 

0  2*1 127E+00 

0  8*222*E*0* 

207  340 

1  000 

1031.94 

0  200939E+00 

0. 70185 3E+0* 

207. 360 

1. 000 

1238  35 

0.  167**9E+00 

0  58*878E+0* 

207  360 

1  000 

1*84  02 

0.  139541ET00 

0. *87398E+0* 

207. 360 

1  000 

1783  22 

0  116284ET00 

0  40616SETO* 

207  360 

1  000 

2139. 86 

0.  969033E-01 

0  338471ET0* 
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TABLE  IS.  PARAMETER  EWH/P  FOR  VARIOUS  R/H  AND  R/L  RATIOS  FOR  A  SIMPLY  SUPPORTED  CYLIN¬ 
DRICAL  SHELL  SUBJECT  TO  A  CENTRALLY  LOCATED  INWARD  POINT  LOAD  BY  FOUR 
ANALYTICAL  METHODS 


E4U4H/P 


r/h  s  a 

1008006+03 

RAOIUS 

THICKNESS 

LENGTH 

iea  oaa 

1.000 

2ae  00 

taa. 000 

1  000 

240  00 

100. 000 

1  000 

288. 00 

100  aea 

i  aaa 

345.40 

iao. aaa 

i.  000 

414. 72 

laa. aaa 

i  000 

497. 44 

100  080 

1  380 

597.20 

100  aaa 

:  aaa 

714. 44 

taa  aaa 

1 . 000 

859  94 

iea. aaa 

1  000 

1031.94 

R/H  3  0 

120000E+03 

120. 000 

1  000 

240  00 

120  oao 

1. 000 

288. 00 

120  000 

1  000 

345.40 

120. 000 

1  000 

414. 72 

120  000 

1  000 

497  44 

120. 000 

1. 00a 

597. 20 

120  aao 

1  000 

714.  44 

120  000 

1. 000 

8S9. 94 

120. 000 

1  000 

1031 . 94 

120.000 

1.000 

1238. 35 

R/H  s  0. 

144O00E+a3 

144  000 

1. 000 

288  00 

144. 000 

1. 000 

345  60 

144.000 

1.000 

414. 72 

144  000 

1  000 

497  64 

144.000 

1.000 

597  20 

144. 000 

1. 000 

714. 44 

144. 000 

1.000 

859.94 

144.000 

1. 000 

1031. 94 

144. 000 

1.000 

1238  35 

144. 000 

1.  BOO 

1484  02 

R/H  s  0 

1728006+03 

172.800 

1. 000 

345 . 60 

172.800 

1  000 

414.  72 

172  800 

1  000 

497  66 

172. 800 

1. 000 

597  20 

172.  800 

1.000 

716.64 

172. 800 

1.000 

859  96 

172.800 

1. 000 

1031  94 

172. 300 

1. 000 

1233  35 

172.800 

1.000 

1486  02 

172.800 

1.000 

1783.22 

R/H  s  a 

2073606+03 

207  360 

1. 000 

414. 72 

207. 360 

1.000 

497  66 

207  340 

1.000 

597. 20 

207  360 

1. 000 

716.64 

£07. 360 

1. 00a 

859  96 

207  360 

1 . 000 

1031  96 

207  360 

1 . 000 

1238. 35 

207  340 

1.000 

1486. 02 

207  340 

1 . 000 

1783. 22 

207  360 

1.000 

2139  96 

R/L 


•. 5000006+00 

a.  4144476+00 

a. 3472226+00 
0. 2B93S2E+O0 
8. 241 1276+00 
a. 200939E+00 
a. 1 674496+00 
a.  1 395  4 16+00 
0  1 16284E+00 
a.  969033E-01 


a. 500000E+00 

a. 416667E+0O 
a. 347222E+00 
a  2893S2E+00 
0  2411276+00 
a. 2009  39E  +  00 
0. 1674496+00 
0. 139541E+00 
0. 116284E+00 
B. 9690336-01 


a. 500000E+00 
a. 414447E+00 
a. 347222E+00 
8. 289352E+00 
8. 241127E+00 
a. 200939E+00 
a.  1 47  4496+00 
a  1  395  4 1  E-f  00 
a. 114284E+00 
a. 949033E-01 


0  500000E+00 
a  4166676+00 
a. 347222E+ea 
a  289352E+00 
a. 241127E+00 
0. 200939E+00 
a  147449E+00 
0. 139541E+00 
0. 114284E+00 
0. 949033E-01 


0. 500000E+00 
0. 4166676+00 
0. 3472226+00 
a. 2893526+00 
a.  2411276+00 
0. 200939E+00 
0. 147449E+00 
0. 1 395416+00 
a. 1162846+00 
8. 949033E-01 


ODGVIST 


0.  143521E+03 
0. 179128E+03 
0. 19422SE+03 
0. 214954E+03 
0. 23S47OE+03 
0  2579446+03 
0. 282564E+03 
0  309533E  +  03 
0. 3390776+03 
a. 371440E+03 


0  205  374E+03 
0  2249786+03 
0  24445 16+03 
0  249974E+03 
0. 295741E+03 
0  323948E+03 
0. 354889E+03 
0  388742E+03 
0  4258676+03 
0. 464514E+03 


0.2579446+03 
0  282564E+03 
0.  3095336+03 
a  339077E+03 
0  371440E+03 
a  404892E+03 
a. 445728E+03 
a  4882706+03 
0. 5348736+03 
0  5B5924E+03 


0  323948E+03 
0  354889E+03 
0  3887626+03 
0  425847E+03 
0  4445 14E+03 
0. 5 1 10416+03 
a  5598176+03 
0. 613249E+03 
0  471781E+03 
8. 735899E+03 


0  404892E+03 
0  445728E+03 
0  488270E+03 
0  5348736+03 
0  S8S9246+03 
0  441848E+03 
0  7031096+03 
a  770218E+03 
0  843731E+03 
a. 924261E+03 


51 JLAARD 


a. 162048E+03 
0.  17B737E+03 
0. 1971756+03 
B. 217584E+03 
a  240179E+03 
a. 245554E+03 
0  294 1 37E+03 
0. 3251796+03 
0. 3S9142E+03 
a  400743E703 


0  202839E+03 
0  2236376+03 
0  244400E+03 
0  2720186+03 
B  3001016+03 
a. 331 309E+03 
0  364543E+03 
0. 405402E+03 
0. 447 471 £303 
0. 49S447E+03 


0. 253794E+03 
0. 2797576+03 
a  306 382E+B3 
0. 340014E+03 
0  374984E+03 
0. 41 35706+0  3 
0  4S47B2E+03 
0  50S387E+03 
a  557971E+83 
0.  615218E+03 


0  317436E+03 
0  349840E+03 
0  385573E+03 
0  424  947E+03 
0  4484S8E+03 
a  5163956+03 
0  549445E+03 
0  4292 34E+03 
0  495  3502+03 
0.7660456+03 


0. 396867E+03 
0  4373726+03 
0  48 1 944E  +  0  3 
0  5  310056+03 
0  585093E+03 
0  444722E+03 
0  710239E+03 
0  7833606+03 
0. 84540 3E+03 
0  9544446+03 


MIZOGUCHI 


0. 141 3 10E  +  0  3 
a. 177950E+03 
a. 194348E+03 
a. 214498E+03 
0. 239194E+03 
a. 24441 3E+83 
0. 292725E+03 
0. 3233S0E+03 
0. 3S4709E+03 
0. 39711 4E+03 


a  2020 44E+03 
0  2228156+03 
0.  24S734E+03 
a  271092E+03 
a. 299074E+03 
0  330125E+03 
B. 345 10SE+03 
0  4037036+03 
0. 444917E+03 
0. 491830E+03 


0. 2S2982E+03 
0. 278917E+03 
0. 307499E+03 
0. 339844E+03 
0. 373933E+03 
0  412357E+03 
0  455 302E+03 
0. 503451E+03 
a  555  326E  +  03 
0. 41 1528E+03 


0. 314434E+03 
0  349027E+03 
0. 384497E+03 
0. 4240046+03 
0. 447409E+03 
0. 515179E+03 
0  S47941E+03 
a. 427298E+03 
0  492442E+03 
0. 742249E+03 


a  3941 17E+03 
0. 4345  90E+03 
0  481 1 17E+03 
0  5  30109E+03 
0  5840876+03 
0  443543E+03 
0. 708788E+03 
0  781443E+03 
0. 8429376+03 
0  950790E+03 


YUAN-TING 


0. 142372E+03 
3. 179509E+03 
0. 199729E+03 
0  221548E+03 
0. 241774E+03 
0.259341E+03 
0  278554E+03 
0.  307784E  +  0  3 
0  3S4590E+03 
0. 4325746+03 


0  20421SE+03 
0  224874E+03 
0  24939SE+03 
0  277448E+03 
0  30S972E+03 
0  3305  39E+0  3 
0. 35  3293E  +  03 
0  383149E+03 
0. 432SS9E+03 
0  S13S89E+03 


0  2S4998E+03 
0. 2823146+03 
0. 311453E+03 
a. 3447056+03 
0  384919E+03 
0  420121E+03 
0. 4500936+03 
0  +82342E+03 
0  531848E+03 
0  4155796+03 


0  323303E+03 
a  3549336+03 
0  390303E+03 
a  4325 196+03 
0. 481702E+03 
0  S31183E+03 
0  S73438E+03 
0  412188E+03 
0  4424316+03 
0. 7470506+03 


0. 404343E+03 
0. 444549E+03 
0  489949E+03 
0.5402576+03 
0  400915E+03 
0  447S42E+03 
0. 7287766+03 
0  7801536+03 
0  3347306+03 

0  919082E+03 
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TABLE  16.  PARAMETER  EWH/P  FOR  VARIOUS  R/H  AND  R/L  RATIOS  FOR  A  PINCHED  CYLINDRICAL 
SHELL  BY  TIMOSHENKO'S  INEXTENSIONAL  SOLUTION 


r/h  s  a 

E*U*H/P 

RADIUS 

THICKNESS 

length 

r/l  iue::tensional 

£90 

see 

l 

000 

S00 

00 

0  1000O0E+01 

0 

324932E+05 

200 

000 

1 

000 

£40 

00 

0. 833333E+00 

0 

270777E+05 

200 

000 

1 

000 

£83 

00 

0  694444E+00 

0. 

225  647E+05 

200 

000 

1 

000 

34S. 

60 

0  S78704E+00 

0 

l  S80  39E  +  05 

200. 

000 

1 

000 

414 

72 

0. 48225 3E+00 

0- 

156699E+05 

200 

000 

1 

000 

497 

66 

0  401878E+00 

0 

130S83E+0S 

200 

000 

1 

000 

S97 

20 

0. 334693E+00 

0 

108S19E+05 

200 

000 

1 

000 

714 

66 

0  279083E+00 

0 

90682684-04 

200 

000 

1 

000 

85  9 

96 

0  £3£S 63E+00 

0 

755638E>04 

200. 

000 

1 

000 

1031 

96 

0  193807E+O0 

0 

629740E+04 

R/H  3  0 

240000E+03 

£40 

000 

1  000 

240 

00 

0  100000E+01 

0 

467902E+05 

£40 

000 

1  000 

238 

00 

0  833333E400 

0 

3899 18E+0S 

£40 

000 

1. 000 

345 

60 

0. 694444E+00 

0 

324932E+05 

£40 

000 

1 . 000 

414 

72 

0  5 78704E  +  00 

0 

270777E+05 

240 

000 

1  000 

497 

66 

0  48225 3E+00 

0 

22S647E4-05 

240 

000 

1  000 

597 

20 

0. 40187SE+00 

0 

138039E4-05 

£40 

000 

1  000 

716 

64 

0.  334898E+00 

0. 

156699E705 

£40 

000 

1. 000 

859 

96 

0. 279082E+00 

0 

130583E+05 

£40 

000 

1.000 

1031 

96 

0. 232S68E4O0 

0 

108819E705 

£40 

000 

1.000 

1238 

35 

0. 193807E+00 

0 

906826E704 

R/H  =  0 

288400E+03 

£88.000 

1 

000 

£38 

00 

0  100000E+01 

0 

673779E+0S 

£88  000 

1 

000 

345 

60 

0  8  3  3  3  3  3E+ 00 

0 

S61483E+0S 

£88  000 

1 

000 

414 

72 

0.  69444 4E 600 

0 

46790EE+05 

£83  000 

1 

000 

497 

66 

0  S78704E+O0 

0 

389918E+0S 

£88  000 

1 

000 

597 

20 

0. 48225 3E+00 

0 

324932E+05 

£83  000 

1 

000 

716. 

64 

0  401878E+00 

0 

27O777E+0S 

£88  000 

1 

000 

8S9 

96 

0  3  34898E+00 

0 

225647E+0S 

£8S  000 

1 

000 

1031 

96 

0.  279082E+00 

0 

188039E+05 

£ca  000 

1. 

000 

1233 

35 

0.  £32St>eE+00 

0 

156699E+05 

£88  000 

1 

000 

1486 

02 

9. 193807ET00 

0 

130S83E+0S 

R/H  =  0 

345  400E+03 

345 

600 

1  000 

345 

60 

0  100000E+01 

0 

970242E+0S 

345 

600 

1  000 

414 

72 

0  833333ET00 

0 

8085  35E  +  05 

345 

600 

1  000 

497 

66 

0  694444E+00 

0 

673779E+0S 

345 

600 

1. 000 

597 

20 

0. S787O4E+00 

0 

561433E+0S 

345 

600 

1  000 

716 

64 

0  48225  3E+00 

0 

467982E+05 

34S 

600 

1  000 

35  9 

96 

0  401878EKO0 

0 

389918E+0S 

345 

600 

1  000 

1031 

96 

0  334898E+00 

0 

324932E+0S 

345 

600 

1.  900 

1238 

35 

•  2790S2EK00 

0 

270777E+05 

345 

600 

1  000 

1486 

02 

0  2  325  oSE  +  00 

0 

22S447E+3S 

345 

600 

1  000 

1783 

22 

0  193807EK00 

0 

188039E+0S 

R/H  =  0 

414720E+03 

414 

720 

1 

000 

414 

72 

0 

1  000  03  £4-0 1 

0 

1 3°7 1SE+06 

414 

720 

000 

497 

66 

0 

8  3  3  3  3  3E+00 

0 

1 16429E+06 

414 

720 

000 

597 

20 

0 

694444E+00 

0 

970242E+0S 

414 

720 

1 

000 

716 

64 

0 

S  78704E+00 

0 

8085  35E  +  0S 

414 

720 

000 

859 

96 

0 

48225  3E+00 

0 

673779E+05 

414 

720 

000 

1051 

96 

0 

401 878E+00 

0 

S61483E40S 

414 

720 

1 

000 

1258 

35 

0 

334898E4-00 

0 

467902E+0S 

414 

720 

1 

000 

1486 

02 

0 

279082E+00 

0 

389918E+05 

*14 

720 

000 

1793 

22 

0 

2325  d8E+00 

0 

324932E+0S 

414 

720 

1 

000 

2139 

86 

0 

193B07E+00 

0 

270777E+05 

46 


NSWC  TR  84-342 


TABLE  17.  PARAMETER  EWH/P  FOR  VARIOUS  R/H  AND  R/L  RATIOS  FOR  A  SIMPLY  SUPPORTED  CYLIN¬ 
DRICAL  SHELL  SUBJECT  TO  A  CENTRALLY  LOCATED  INWARD  POINT  LOAD  BY  FOUR 
ANALYTICAL  METHODS 

E6U6H/P 


R/H  =  0  200000E+03 


RADIUS 

THICKNESS 

LENGTH 

R/L 

0D0VI5T 

B  t  JLAARD 

HI ZCGUCHI 

200  000 

1. 000 

200  00 

0. 100000E+01 

0. 27S003E  +  03 

0. 24163AE+03 

0. 260958E+03 

280  000 

1  080 

240.00 

0. 833333E400 

0. 301256E+03 

0  2887S2E+03 

0  2680S0E+03 

200. 000 

1. 000 

288. 00 

0  694444E  f 00 

0. 330010E+03 

0.  318SO5E  +  03 

0. 317781E+03 

200  000 

1  000 

345 . 60 

0. 578704E+00 

0.  361507E+03 

0. 3S1179E+03 

0.  35  04  33E+03 

200. 000 

1  000 

414. 72 

0.  48225 3E+00 

0.  396012E+83 

8  387083E+03 

0  386  3 13E  +  03 

200  000 

1  000 

497. 66 

0.  401S78E+00 

0. 433809E403 

0. 42653SE+03 

0. 425732E+03 

200.000 

1.000 

597  20 

0. 33439SE+00 

0. 475214E+03 

0. 4700O1E+03 

0.  469209E+03 

200  000 

1 . 000 

716  64 

0. 279082E+00 

0  520571E+03 

0  S17937E+03 

0. 5 17010E+03 

200  000 

1. 000 

859  96 

0. 232S68E+00 

0. S70257E+03 

0.  57O770E+03 

0. 569724E+03 

200.000 

1. 000 

1031.96 

0. 193807E+0O 

0. 62468SE+03 

0. 628952E+03 

0. 6277 13E+03 

R/H  a  0 

2»0O00E+O3 

240  000 

1 . 000 

240  00 

0. 100000E701 

0  345400E+03 

0  326771E+03 

0. 3261 94E+02 

240  000 

1 . 000 

288  00 

0. 8  3  3  3  3  3E  +  00 

0  37836&E+03 

0.  360724E403 

0. 360124E+03 

240  000 

1 . 000 

34S . 60 

0  694444E+0© 

0. 414480E+03 

0.  397968E+03 

0. 397346E+02 

240  000 

1  000 

414. 72 

0. 5  7S704E+00 

0. 4S4040E+03 

0.  438849E  +  83 

0. 436206E+03 

240  000 

1 . 000 

497. 66 

0. 48225  3E+00 

0. 497375E+03 

0  483752E+03 

0  48308  3E  +  0  3 

240 . 000 

1. 000 

597.  20 

0. 40 1 978E+00 

0  5 44848E+0 3 

0  S  33098E+03 

0  S32395E+03 

240. 000 

1. 000 

716. 64 

0  334S98E+00 

0  59665 1E+03 

0  S87357E+03 

0  586603E+03 

240. 000 

1 . 000 

859  96 

0. 279082E+00 

0  65  3817E  +  03 

0  647022E403 

0. 646 189E+03 

240  000 

1 . 000 

1031.96 

0. 2325  &8E+00 

0. 71 6221E+03 

0. 712679E+03 

0  71 1723E+03 

240  000 

1  000 

1238. 35 

0.  193807E+00 

0  7845 80E  +  03 

0. 785  80S  E  +  0  3 

0  78 36S 5 E  +  0 3 

R/H  3  0. 

288000E+03 

283.000 

1.000 

283.  00 

0.  10OO0OE+01 

0.  433809E403 

0. 407722E+03 

0  407336E+O3 

283.000 

1. 000 

345  60 

0  833333E+00 

0. 475214E+03 

0  4S0246E403 

0  449837E+03 

283.000 

1.000 

414.72 

0  694444E+00 

0. S20571E+03 

0.  496873E+03 

0  496449E+03 

233. 000 

1  000 

497  66 

0. S78704E+00 

0. 57025 7E403 

0  S48043E+03 

0  547S91E+03 

288  000 

1  000 

597  20 

0  48225 3E+00 

0. 624685E+03 

0  604204E+03 

0  4O3728E+03 

233  000 

1.000 

716  64 

0. 401878E+00 

0  634308E+03 

0  665377E+03 

0  665  362E+03 

288.000 

1  000 

85  9  96 

0. 334898E+00 

0. 749622KT03 

0  733607E403 

0  733038E+03 

238  000 

1 . 000 

1031  96 

0. 279082ET00 

0  821 1 69E  +  0 3 

0.  8O7993E+03 

0  BO7342E+03 

238  000 

1  000 

1233  35 

0. 232S68E+O0 

0  399546E+03 

0  3S964HE+03 

0  383S69E+03 

283.000 

1  000 

I486  02 

0. 193807E+00 

0. 985403E+03 

0. 979429E+03 

0. 973452ET03 

R/H  =  0 

345600E+02 

345 

600 

1  . 

000 

345 

60 

0  100000E+01 

0 

544848E+03 

0  50815 9E+03 

0 

5  0S09SE+03 

345 

600 

1 

000 

414. 

72 

0. S33333E+00 

0 

596851E+03 

0  561429E+03 

0 

S61345E+03 

345 

600 

1  . 

000 

497 

66 

0. 694444E+00 

0 

65  3817E+03 

0. 619829E+03 

0 

6 1 9723E  +  0  3 

345 

.  600 

1 

000 

597 

20 

0. 578704E+00 

0. 

716221E+03 

0. 68387SE+03 

0 

633749E+0  3 

345 

600 

1 

000 

716 

64 

0. 48225  3E+00 

0 

784580E+03 

0. 75  4141 E+03 

0 

75  39Q4E  +  03 

345 

600 

1 

000 

85  9  . 

96 

0. 401378E+00 

0 

8S9465E+03 

0  831 232E+02 

0 

S31036E+03 

345 

600 

1 

000 

1031 

96 

0.  3  34S98E+00 

0 

941497E+03 

0  9 1£ 800E+03 

0 

91  SS47E+03 

345 

.  600 

1 

000 

1238 

35 

0  279082E+00 

0. 

103136E+O4 

0  10*385  5  E  +  04 

0 

100S21E  +  04. 

345 

.  600 

1  . 

000 

1486 

02 

0. 2325  68E+00 

0 

1 1 29  8  0E+04 

0  11101 9E+ 04 

0 

1  10971E+04 

345 

.  600 

1. 

000 

1783. 

22 

0. 192807ET00 

0 

123763E+04 

0  122161E+04 

0. 

122093E+04 

R/l 

K  =  0 

414720E+03 

414 

720 

1 . 000 

414 

72 

0. 1 00000E+0 1 

0 

A84308E+03 

0 

6 3ri5  32E+03 

0 

6  32988E  +  0 3 

414 

720 

1 . 000 

497 

66 

0. 8  33  3  3  3E+00 

0 

749422E+03 

0 

699277E  +  0  3 

0 

699709E+03 

414 

720 

1 . 000 

5  97 

20 

0  694444E400 

0 

92U69ET03 

0. 

772427Et03 

0 

7723  36E+03 

414 

720 

1  000 

716 

64 

0. S78704E+00 

0 

8995A4E+03 

0 

852618E4-03 

0 

8S3O04E+03 

41i 

720 

i  00 : 

359 

96 

0. 4822S3E+00 

0 

985A03E+03 

0 

4405  3 ?E+03 

0 

940894E403 

414 

720 

l.  000 

1031 

96 

0. 401878E+00 

0 

107946E+0* 

0 

10  3692E+04 

0 

103723E+Q4 

414 

720 

1 . 000 

1233 

35 

0  334898E400 

0 

1 1S243E+Q4 

0. 

114254E404 

0. 

114274E404 

414. 

720 

1  000 

1486 

02 

0  279082E4O0 

0 

1295  35E+04 

0 

125319E+04 

0 

1 25  835  E+04 

414 

720 

1  000 

1733 

22 

O. 2  325  63E+0O 

0 

1A1398E+0* 

0 

138476E+04 

0 

138477e+®4 

414 

720 

1  000 

2139. 

36 

0  193807E+00 

0. 

15S4*aE+0* 

0 

152319E+04 

0 

152294E404 

YUAN-TING 


0  269812E+03 
0. 296931E+03 
0  326685E+03 
0  359593E+03 
0  395784E+03 
0. 434608E+03 
0. 47720&E+03 
0. S27674E+03 
0. S87668E+03 
0  6S0  35  6E+03 


0. 339S21E+03 
0. 3735 15E+03 
0  41075 4E+03 
0  4S 1 765E+03 
0  49752&E+03 
0  S46S92E+03 
0  5994oOE+03 
0. 6S9703E+03 
0  732619E+03 
0  SIS  3 1SE+03 


0  427143E+03 
0  469734E+03 
0  516422E+03 
0. S67636E+03 
0  624344E+03 
0  686876E4-03 
0.  75  3832E+-03 
0. 826941E+03 
0  912756E+03 
0. 101772E4O4 


0  5372S0E+03 
0  590627E+02 
0. 649133E+02 
0. 7 1 333 1E+03 
0  783952E+03 
0  S6241^E+03 
0  9477 2EE+03 
0  103877E+04 
0. 1 1 4c2cE+04 
0  1268 10E+04 


0  675  695E+0 3 
0  742533E403 
0  915 903ET03 
0  896398E+03 
0  93461SE+03 
0  108238E+04 
0  119044E+04 
0  1  iQ60lE  +  v54 
0  143173E+04 
0  153122ET04 
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NOMENCLATURE 


Cjj, ,  c|{j  =  Fourier  coefficient 


Eh' 


=  Flexural  rigidity  (=  0 

12(l-v  ) 


) 


E 

h,H* 

J 

L,  J, 


*  Young's  Modulus 

*  Thickness  of  shell 

*  v^r 

*  Length  of  shell 


*  Bending  moraent/unit  length  of  a  section  perpendicular  to  the  axis  of 
the  cylindrical  shell,  i.e.,  per  unit  of  circumference 


M 


X<P 

Mcp 


m 


m ,  n 


=*  Torsional  moment/unit  length  of  an  axial  section  of  cylindrical  shell 
=  Bending  moment/unit  length  of  axial  section 
=  Axial  harmonic  index 
=  Peripheral  harmonic  index 
=  Concentrated  force  applied  radially 
*  Fourier  coefficient 
=  Mean  shell  radius 


RU,RV,RW  =  Rotations  about  local  (U,V,W)  axes  (Figure  1) 


=  Peripheral  coordinate  along  a  circular  arc  on  cylindrical  shell 
(=  Rep) 


=  Axial  coordinate  (Figure  1) 


(U,V,W)  -  Displacements  along  local  axial,  local  tangential,  and  local 

radial  axes  (Figure  1) 


*Except  in  Equation  (12) 
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APPENDIX  A 
SOLUTION  METHODS 


In  this  Appendix,  we  focus  our  attention  to  closed  form  solutions,  which 
are  primarily  based  on  linear  shell  theories, A-^  ~  A-10  for  t]ie  response  of 
cylindrical  shells  subjected  to  static  point  loads.  Small  displacements, 
strains,  rotations,  and  linear  elastic  material  are  assumed  throughout  this 
work.  Solutions  based  on  three-dimensional  elasticity  theory  are  substantially 
more  complex  and  will  only  be  mentioned. 

The  linear  shell  equations  can  be  written  as  a  partial  differential 
equation  of  eighth  order  in  the  radial  displacement  w  by  eliminating  the  other 
two  deformation  measures,  u  and  v.  Such  are  the  full  Fluegge ,  A~*-  simplified 
Donnell, A-^  improved  Donnell ,  A-12  Morley ,  A_l-* »  A-^  etc.,  shell 
equations.  If  certain  terras  are  neglected,  the  range  of  applicability  can  be 
impaired.  In  Appendix  C  it  is  demonstrated  how  Timoshenko  equationsA_2  lead 
to  a  nonsymmetric  coefficient  matrix  in  contrast  to  those  of  VlasovA-15  or 
Fluegge . 

Most  of  the  solutions  obtained  in  the  literature  are  based  on  the 
simplified  DounellA-H-»  A~^  equation,  rather  than  the  complete  Donnell, A-^^ 
the  full  Fluegge ,A_^  Morley,A-“«  A”^  Schorer  ,A~^  >  A-18  etc., 
equations .  ^There  are  numerous  solutions  based  on  Novozhilov ,A"9 

and  other  theories  critically  analyzed  in  Buchwald . A_^9 ,A~20 
20  it  is  shown,  through  asymptotic  expansions  that  different 
solutions  in  the  radial  displacement  w  are  required  for  an  "inner"  and  "outer" 
regions.  Some  of  the  simplified  equations  such  as  Donnell ' sA-H »  A-16  or 
Morley ' sA_^ *  A-^  attempt  to  cover  the  entire  range.  They  contain  incorrect 
second  order  terms  in  both  regions,  and,  therefore,  in  limiting  cases  the 
dominating  behavior  may  not  be  correct.  In  each  region  the  second  order  terms 
must  not  be  included  in  the  approximate  theories  for  limiting  behavior,  although 
the  first  order  solution  over  the  entire  cylinder  may  be  correct. 

Before  reviewing  the  pertinent  literature,  a  few  words  are  needed  to 
establish  the  general  framework  of  solution  methods  presently  available  to 
cylindrical  shells  (Figure  1  of  the  main  text)  subject  to  a  variety  of  loads  and 
end  boundary  conditions.  The  theory  of  variational  or  energy  based  techniques 
(such  as  finite  element  method  (FEM))A-21  will  not  be  referred  to  here. 

The  analytical  methods  based  on  shell  theory  can  be  grouped  in  the 
following  classes,  irrespective  of  the  partial  differential  equation,  which  must 
be  solved. 


Sanders,  . 
ThereA_l9*  A~ 


A- 1 
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SEPARATION  OF  VARIABLES 

Essentially  most  of  the  shell  solutions  are  of  this  form,  except  for  the 
ones  by  asymptotic  methods.  The  method  of  solution  may  proceed  under  either  1, 
2,  3,  or  4  below: 


1.  Resolve  the  loading  function  into  a  double  Fourier  series  of  the  form 

Pmnsin(mirx/L)  sin(n<p) 

ra  n 


for  simply  supported  ends.^-^*  A-23  xhe  choice  of  trigonometric  products 
must  be  such  that  the  boundary  conditions  are  satisfied  automatically.  This 
method  was  originally  given  by  H.  Reissner^-^  for  a  simply  supported 
cylinder.  To  accelerate  convergence  one  may  resort  to  tricks  such  as 
subtracting  and  adding  a  series  corresponding  to  a  reference  plate 
problem. The  solution  of  the  plate  problem  is  obtained  by  other 
techniques.  The  direct  approach  is  to  incorporate  a  large  number  of  terms  in 
the  series.  If  this  last  route  is  followed,  calculation  of  stress  resultants 
and  couples  becomes  less  accurate  than  displacements  due  to  the  term  by  term 
successive  differentiations. 


2.  Resolve  the  loading  function  and  displacements  in  a  single  Fourier 
series  peripherally  and  exponential  functions  axially  of  the  form 


exp(Xc ) 


sin(m<p) . 


The  representation  by  a  Fourier  series  peripherally  and  Fourier  integral  axially 
is  grouped  under  this  class  too.  It  is  worth  noting  that  slowly  damped  edge 
disturbances  must  be  included  to  complete  the  solution.  This  method  is  suitable 
to  treat  line  loads  on  generatrices  or  loads  applied  at  edges  x=const. 

3.  Resolve  the  radial  displacement  w  in  a  single  Fourier  series  axially 
and  damped  trigonometric  function  peripherally  of  the  form 


sin(.-j— )  or 


£ 


,C  exp(n<p)  sin(n<p)  sin(^~) 

1  n  L 


Edge  disturbances  from  the  generatrix  are  needed  to  complete  the  solution.  This 
method  is  better  suited  to  treating  loads  at  the  edges  <p=const. 

4.  Obtain  the  particular  integral,  oj(x,<p),  using  1  above.  Resolve 
the  solution  of  the  homogeneous  8th  order  equation  in  a  series  of  the  form 


00 

PTT(x,<p)  = 


P  (x)  cos  (m(q>  +  a  )  ] 
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where  pm(x)  is  the  solution  of  the  8th  order  linear  homogenous  ordinary 
differential  equation  with  constant  coefficients  and  satisfies  the  boundary 
conditions ,A“24  similar  stress  function  methods  can  be  found  in  References 
A-15,  A-25,  and  A-26. 

As  an  additional  comment,  we  note  that,  if  a  solution  to  Donnell's  shell 
equations  is  known,  then  we  may  obtain  solutions  of  the  full  Fluegge  equation  by 
using  1,  2,  3,  or  4.  Correction  functions A“27  are  obtained  by  the  previous 
method  of  separation  of  variables. 


TRANSFORM  METHODS  (A  form  of  Separation  of  Variables) 

If  the  equations  of  elasticity  are  employed,  then  the  infinitely  long  shell 
can  be  treated  by  the  following: 

a.  Fourier  integral  axially,  Fourier  series  peripherally,  and  stress 
functions  methods ,A~28 

b.  Fourier  integral  transform  methods  with  associated  stress 
function. A-29 


ASYMPTOTIC  INTEGRATION  METHODS 

These  techniques  can  be  found  in  works  by  E.  ReissnerA-30  and  KoiterA“31 
and  references  listed  therein. 


At  this  stage  it  must  be  mentioned  that  the  eighth  order  algebraic 
characteristic  equation  has  eight  complex  roots.  Of  these,  two  are  the 
so-called  principal  roots.  They  are  such  that  their  real  parts  are  negative  and 
imaginary  parts  are  positive  (second  quadrant  on  complex  plane).  These  are  the 
highly  and  lightly  damped  roots.  However,  for  the  two  particular  values  of  ra 
(m=0  or  m*l)  the  characteristic  equation  has  zero  roots  of  multiplicity  4.  The 
assumed  form 


w(t  ,«p) 


oo 

z 

mss  n 


♦  (t)  cos(m<p)  ■ 
m 


^  ^  'y  *  C^  expCX^^t)  cos(mip) 
m=Q,l  k= 1 


does  not  account  for  all  of  the  zero  roots  of  multiplicity  4.  There  are  12 
elementary  states  of  stress  that  correspond  to  the  zero  roots. A  ^  Six  are 
rigid  body  motions,  while  the  other  six  are  elementary  states  of  stress  with  the 
shell  looked  upon  as  a  beam. 
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HISTORICAL  REVIEW 

Linear  shell  theories  have  been  developed  by  Novozhilov,A_9 
Gol'denveizer ,A-3  E.  Reissner , A-33 »  A-32  Fluegge , A-^  Love, A" 7 
Timoshenko, A~2  Sanders  Jr.,A_1°  Koiter,A'31»  A_33  and  Vlasov. A_l5>  A~25 

The  earlier  work  on  specific  problems  on  circular  cylindrical  shells 
was  carried  out  by  H.  Reissner ,A-22  Miesel , A“34  Fluegge, 
Aas-Jakobsen,A_33»  A-33  Finsterwalder , A_37  Schorer,A~^7  Odqvist , A_^® 
Donnell, A“H  Bi j laard , A-23  Yuan,A_33  Yuan  &  Ting,A_24»  A-39 
Morley,^3^"1^  Hqff,A_40  “  A-^2  Kempner ,A_43 * A“44  Mizoguchi  A“45 
Seide,  ®  Yao,A  23  Bieger,A  47  Berger, A~43  Melworm,  et  al.A-4^  and 
others . 


More  specifically,  FinsterwalderA-37  and  SchorerA_^-7  simplified  the 
various  equations  under  restrictive  assumptions.  SchorerA-^7  reduced  the 
eighth  order  homogeneous  equation  to  a  fourth  order  one  by  neglecting  lower 
order  derivatives.  He  based  this  on  the  assumption  that  damping  factors  exceed 
unity;  therefore,  the  higher  the  derivative,  the  more  dominating  it  becomes. 
This,  of  course,  is  not  always  the  case.  He  also  introduced  additional 
corrections  based  on  Newton's  method  to  improve  upon  his  solution. 

OdqvistA-^3  used  Schorer's  theor^-^7  to  obtain  the  only  truly  explicit 
solution  for  a  simply  supported  circularly  cylindrical  shell  subjected  to  a 
concentrated  load. 


BijlaardA~23  solved  the  circularly  cylindrical  shell  problem  subject  to  a 
variety  of  loads  and  simply  supported  conditions  by  using  a  double  Fourier 
series  method.  The  eighth  order  differential  equation  is  slightly  different 
from  Donnell' sA-^l  due  to  the  simplifying  assumptions.  Limitations  of  the 


analysis  occur  for  certain  combinations  of  L/R  and  R/h,  as  pointed  out  by 
Murphy . A" 2 3 


Mizoguchi ,A-A5  employed  his  own  version  of  the  characteristic  equation. 

He  solved  the  simply  supported  shell  problem  by  a  finite  Fourier  transform 
method  and  the  method  of  images  (Figure  2  of  the  main  text).  YuanA  33  used 
the  simplified  Donnell  equation4-^  to  obtain  a  solution  of  the  infinitely 
long  cylinder  subject  to  two  self-equilibrating  point  loads.  The  load  was 
resolved  in  a  product  of  a  Fourier  integral  axially  and  a  cosine  Fourier  series 
peripherally,  since  the  function  is  even  in  both  longitudinal  and 
circumferential  variables  (x,s).  Ting  &  YuanA  3^  employed  the  complete 
Donnell  equation. A  In  their  analysis  they  treated  boundary  conditions 
other  than  simply  supported.  The  infinite  length  cylinder  solution  was  used  as 
a  particular  integral.  Boundary  conditions  were  satisfied  by  obtaining  the 
proper  solutions  to  the  homogeneous  differential  equation.  Yuan  &  TingA 
applied  Fluegge'sA  full  equations.  Both  loading  function  ana  radial 
displacement  were  resolved  in  a  product  of  a  Fourier  integral  axially  and  a 
Fourier  series  peripherally.  The  displacement  of  the  infinitely  long  cylinder 
is  obtained  first.  They  used  the  method  of  images  to  obtain  the  correct  end 
boundary  conditions  of  a  simply  supported  cylinder  of  finite  length. 
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SeideA“^6  analyzed  the  bending  of  an  infinitely  long  cylindrical  shell 
subject  to  point  loads  or  moments  using  Donnell's  equations.  Similar  work  was 
done  by  Melworm,  et  al.A-^9  BiegerA-^  reduced  the  three  coupled  partial 
differential  equations  of  FlueggeA_^  for  the  three  deformation  measures 
(u,v,w)  to  an  eighth  order  equation  for  an  auxiliary  function  F(x,<p)  for  the 
case  of  a  concentrated  load.  He  solved  the  infinitely  long  cylinder  problem  by 
resolving  the  load  in  a  product  of  a  Fourier  series  peripherally  and  a  Fourier 
integral  axially. 

CooperA-^0  employed  Naghdi's  shallow  shell  equations,  including  the 
effect  of  transverse  shear  deformation,  in  order  to  obtain  solutions  to  the 
simply  supported  cylindrical  shell  of  finite  length  subject  to  a  line  load  on 
the  generatrix.  His  results  indicated  close  agreement  stresswise  with 
DonnellA~H  equation  solutions. 

Biezeno  and  KochA“-^  gave  particular  solutions  to  a  system  of  three 
coupled  partial  differential  equations  for  an  assumed  load  system  radially, 
tangentially,  and  axially.  Boundary  conditions  must  be  accounted  for 
separately.  This  can  also  be  found  in  Reference  A-5. 

BergerA-^®  studied  the  roots  of  the  homogeneous  FlueggeA-l  equation, 
assuming  that  the  displacement  varied  as  expC-Xx/a)  exp(jn*p).  He  gave 
comparisons  with  the  solutions  by  MieselA~-^  and  Aas-Jakobsen. A~35 

This  brings  us  to  the  topic  of  eigenvalues  of  the  characteristic  equation. 
Morley^-  calculated  the  eigenvalues  of  the  homogeneous  differential  equation 
by  Fluegge,A-l  Donne 11,A-^  and  his  own  version.  He  showed  the  proximity  of 
his  to  Fluegge's  for  various  values  of  the  harmonic  index  n  in  contrast  to 
Donnell 1 s ,A-I1  which  deviates  more,  the  thicker  the  shell  and  the  smaller  n 
(either  peripheral  or  axial  harmonic  content)  is.  Furthermore,  special  methods 
were  used  to  solve  the  eighth  order  characteristic  equation  by  reducing  it  to  a 
biquadratic.  For  a  general  biquadratic  there  exist  generalA~->2  and  special 
methods .A-5 3 

MorleyA-^  also  solved  the  inhomogeneous  problem  for  his  own  equation.  A 
similar  study  for  Donnell's  equationA-^^  with  axial  prestress  (tensile  or 
compressive)  was  conducted  by  Nachbar,A~^  whose  analysis  was  very  sketchy 
with  missing  intermediate  results.  Some  can  be  recovered  from  Holand.A~^ 
Furthermore,  there  is  a  series  of  articles  by  HoffA-i!*^>  A “ ^  with  explicit 
solutions  for  Che  roots  of  the  homogeneous  DonnellA-^  differential  equation 
as  well  as  the  membrane  stresses  and  displacements.  It  was  concluded  in 
Reference  A-42  that  under  certain  conditions  Donnell's  equationsA-H  can  be 
used  and  still  be  accurate  in  preference  to  the  more  complicated  FlueggeA-^ 
equations.  Further  discussion  appears  in  Reference  A-44.  However,  the  most 
extensive  study  seems  to  be  the  books  by  HolandA“^  and  Seide,A-^^  where 
other  approximations  to  the  large  and  small  roots,  apart  from  previous  ones, 
such  as  Gruber '  s  ,A~56  Miesel's, Aas-Jakobsen '  s  ,A-35  Berger  '  s  ,A_/f® 

Olsen's, A-^  etc.,  are  given  and  plotted  in  an  attempt  to  show  that,  at  least 
for  some  engineering  applications,  special  savings  can  be  obtained  if  one 
resorts  to  them.  Other  studies  to  be  mentioned  are  the  works  of  Houghton,  et 
al.A“57»  an<j  HolandA-^  and  a  recent  study  on  the  form  of  thin  shell 
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equations  and  their  eigenvalues  by  Akeju.A-60  Furthermore,  analysis  by 
Buchwald^'l^ ,  A-20  indicates  the  importance  of  obtaining  the  correct  exact 
eigenvalues  for  limiting  cases. 

In  view  of  experimental  evidenceA-kl  it  is  not  advisable  to  neglect  the 

supported  cylindrical  shell  problem  subject  to  numerous  loading  conditions  on  a 
generatrix.  The  radial  displacement  w  was  assumed  as  a  Fourier  series  axially 
and  an  exponential  peripherally  with  the  real  part  negative. 

Hoiand,A“62  employing  Donnell ' sA"H  theory  for  a  simply  supported 
cylindrical  shell,  resolved  forces  and  loads  in  a  single  Fourier  series 
longitudinally  and  a  damped  trigonometric  function  peripherally.  He  was 
principally  interested  in  obtaining  influence  surfaces  for  bending  moments. 
Stafford  and  Rira^-^  introduced  a  new  approach  by  a  corrector  function  to  the 
Donnell'^-H  equation  solution  through  an  iterative  procedure. 

BuchwaldA-19 ,A-20  wag  concerned  with  the  various  differential  equations 
derived  from  different  shell  theories  as  applied  to  thin  circular  cylinders. 
There  are  regions  of  validity  for  each  of  them  (Donnell,  Fluegge,  complete 
Donnell,  Morley,  etc.),  i.e.,  inner  and  outer  regions  as  well  as  regions  of 
transition. 

Kempner,  et  al.,A-^^  employed  DonnellA“H  equations  to  obtain  tables 
and  curves  for  simply  supported  shells  subject  to  numerous  loading  conditions 
for  a  variety  of  radius/thickness  and  radius/length  ratios. 

YajA-^  treated  the  problem  of  a  thick  long  cylinder  subject  to  two 
radial  self-equilibrating  loads  by  elasticity  theory.  The  load  is  represented 
by  a  Foirier  series  circumferentially  and  a  Fourier  integral  axially.  A  stress 
function  is  constructed  and  appropriate  boundary  conditions  are  satisfied  (all 
stresses  must  vanish  at  +  ®  &  -<*>)  while  stresses  at  inner  and  outer  radii 
are  prescribed.  The  radial  displacement  is  finally  given  in  terms  of  modified 
Bessel  functions. 

Klosner^-^  solved  the  long  circular  cylindrical  shell  problem  subjected 
to  uniform  external  circumferential,  radial  line  loads  by  employing  elasticity 
equations  and  Fourier  integral  transform  techniques.  The  inverse  transforms 
were  obtained  by  numerical  and  asymptotic  integration. 

Finally  the  solution  of  an  infinitely  long  cylinder,  reinforced  by  a  single 
ring  and  subjected  to  self-equilibrating  point  loads  eccentrically  applied  was 
solved  by  Allentuch,  et  al.A_^,  by  a  Fourier  ser ies-Fourier  integral  method. 
Their  results  were  supported  by  test  data.A-^  The  unusual  point  was  that 
their  radial  deformation  was  not  a  maximum  under  the  ring  stiffener.  It  was 
claimed  to  be  the  effect  of  eccentricity.  To  this  we  must  add  the  inextens iona l 
solution  by  TimoshenkoA-^  (sometimes  referred  to  as  Rayleigh's  solution)  of  a 
circularly  cylindrical  shell  subject  to  point  loads.  Note  that  the  ends  are 
free  and  there  is  no  membrane  action  in  the  longitudinal  direction  of  the  shell. 
This  solution  resembles  the  ring  subject  to  point  loads  solution ,A-65 


longitudinal  moment,  Mjj,  and  torsional  moment,  M,™, .  Consequently,  Hoff 
et  al.^-1^  employed  Donnell  equationsA“H »  A-l6  to  soive  the  simply 
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APPENDIX  B 

BIJLAARD'S  SOLUTION  BY  DOUBLE  FOURIER  SERIES 

CIRCULAR  CYLINDER  OF  FINITE  LENGTH  SUBJECT  TO  A  POINT 
LOAD  AT  MIDSPAN  SIMPLY  SUPPORTED  AT  BOTH  ENDS 


The  following  presentation  is  based  on  Bijlaard's  analysis. We  employ 
Timoshenko's  equations**"^  (eq.  303,  p.  513)  of  equilibrium  for  the  three 
deformation  measures  (u,v,w),  allowing  only  radial  pressure  load,  pr,  with  E, 
Young's  Modulus,  v,  Poisson's  ratio,  R,  the  cylinder's  mean  radius,  h,  the 
skin  thickness  (Figure  B-l  shows  a  point  load  P  applied  at  the  cylindrical 
shell) . 
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where  x  is  the  longitudinal  shell  axis,  tp  the  circumferential  direction,  and  r 
the  radial  axis.  The  object  of  the  following  paragraphs  is  to  eliminate  the 
displacements  u,v  and  obtain  a  partial  differential  equation  for  the  radial 
deformation  w,  only  involving  the  loading  function  pr  also.  Because  of  the 
process  of  repeated  differentiation,  not  all  of  the  solutions  of  the  reduced 
eighth  order  equation  satisfy  equations  (B.l),  (B.2),  and  (B.3).  Additional 
solutions  enter,  which  must  be  excluded  by  careful  analysis. 

The  final  partial  differential  equation,  which  was  obtained  after  a  lot  of 
algebra  and  neglection  of  terms  of  order  [h/R]^  or  higher,  assumes  the  form 


„8  A  12(l-v  ) 
V  w  + 


rV 
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a  w  .  i_ 
ax4  r2 


3  ^w  +  (\>  +  7 )  3  ^w 
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36w 
3  x43  <p  2 
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D  rr 
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Eh' 


12(l-v  ) 


Equation  (B.5)  was  solved  by  the  method  of  separation  of  variables  through 
Fourier  series  (the  fact  that  the  P.D.E.  contains  even  powers  of  w  derivatives 
is  decisive).  In  particular  we  may  represent  the  radial  load  by  a  double 
Fourier  series  and  similarly  the  radial  displacement.  Such  an  operation  leads 
to  a  term  by  term  equality  of  the  coefficient  of  the  series.  Note,  however, 
that,  by  resolving  the  load  in  such  a  form,  boundary  conditions  and  any  symmetry 
or  antisymmetry  considerations  must  be  satisfied.  Unlike  ordinary  differential 
equations,  where  the  homogeneous  (complementary  function)  and  particular 
solutions  (particular  integral)  produce  the  solution,  partial  differential 
equations  are  substantially  more  difficult  to  solve  and  obtain  general 
solutions.  We  are  content  to  find  just  a  solution. 


NSWC  TR  84-342 


Let  the  surface  loading  function  pr  be  resolved  in  a  series  of  the  form 
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where 
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To  obtain  the  Fourier  coefficients  we  may  proceed  as  follows.  By  separation 
of  variables  resolve  pr  (x,s)  into  a  product 
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q(s )  =  Qo 

where 

q(s)  =  0 


f°r  |  s  |  <_  Cj. 
otherwise 


(B.16; 


p(x)  =  P0 

and 

p(x)  =  0 

and  b  m  L/2,  where  L  is  the  entire  length  of  the  shell.  Since  q(s)  is  even  in  the 
arc  variable  s,  we  may  resolve  it  in  a  cosine  Fourier  series  in  the  interval 
(-ira,  na)  with  a  period  T  =  2ira,  i.e., 
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Similarly,  p(x)  is  resolved  in  an  odd  Fourier  series  in  the  interval  (-L,+L)  with 
period  2L  and  b=L/2 


p(x)  =  £  Bn 

n=l 


(B.22] 
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Finally,  the  point  load  case  can  be  arrived  at  by  taking  lirats  as  both  and  6 
tend  to  zero  and 


4Q0P0ci6  =  P 
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The  Fourier  coefficients  Zm  n  for  the  loading  function  of  a  radially 
concentrated  load  at  the  midpoint  of  span  L  (even  series  peripherally)  are 
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Consequently,  the  radial  deflection  finally  can  be  written  as 


w  w 
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MIZOGUCHI's  SOLUTION 


Mizoguchi®-3  employed  a  slightly  different  form  of  the  eighth  order 
partial  differential  equation  in  radial  displacement  w  (again  even  powers  of  w 
derivatives  of  w) 
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He  applied  a  double  Fourier  transform  to  both  sides  of  (B.30)  of  the  form 

w(m,n)  =  J  J  w(x,0)  cos  — r-  cos  (nq>)  d<pdx  (B.3 

-L  -it 

He  considered  the  cylindrical  shell  to  be  subjected  to  a  series  of  point  loads 
(Figure  B-2),  alternating  from  +P  to  -P  every  L.  Consequently  at  +L./2 

w=0,  v=0,  Nx  =0,^=0  (b.3 

At  jHL,  where  the  deflection  is  maximum,  the  slope  is  zero.  To  use  this  method 

a3  a5 

o  W  o  W 

we  need  to  know  the  values  of  the  boundary  terms  — —  and  — —  at  +L/2. 
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On  the  other  hand, 

p(x,<p)  *  £  y. 
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if  we  proceed  as  Bijlaard 
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Therefore,  the  displacement  under  the  point  load  P  becomes,  noting  that,  when 
the  peripheral  index  n=0,  the  numerator  and  denominator  simplifies 


(B. 38) 


n-1 

2 

where  (-1)  sin(nni/2)  has  been  replaced  by  1.0  for  the  odd  values  that  ra  and 
n  assume. 


B-9 


NSWC  TR  84-342 


ODQVIST's  SOLUTION 


Odqvist®-^  used  Schorer's  theory®--^  to  obtain  the  solution  to  a  simply 
supported  shell  subject  to  a  concentrated  load  P  at  midspan  (also  line  load  on 
generatrix).  Starting  from  the  simplified  equation 


12(1-v2)R2 


(] 


and  assuming  a  general  term  of  the  Fourier  series  to  be 


he  determined  in  the  next  steps  the  form  of  the  series.  He  used  the  principal 
roots  (second  quadrant  as  mentioned  previously)  from  the  auxiliary  equation 
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Finally , 


the  form  of  the  Bending  Moment/unit  length 


is 


(B 


Similar  assumptions  for  the  load  and  a  lot  of  manipulation  yields  that  the 
central  radial  displacement  is 
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Because  the  computation  of  this  infinite  series  is  extremely  slow  (over  20,000 
terms)  we  may  proceed  as  follows  to  evaluate  the  series. 

From  Whittaker  &  Watson®-®  we  have  that  the  definition  of  the  Zeta 
function  is 


c(s)  =  l/nS  ( 

n-1 } 2  j  t • * 

where  s  could  be  complex.  This  is  an  analytic,  uniformly  convergent  series  for 
Re(s)  >1+6 .  Furthermore,  by  the  Jensen  formula,®-®  (p.  279) 
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We  note  that  as  6*0,  the  integrand  -*-s 6/2,  while  the  portion  of  the 
integrand 


[1  +  exp  (ir  tan  8)]  cos^e  + 


ir 


2 


2 


However , 


-  -4  as) 

2  n=l 


■L 


i 

(2n)S 


y  — L— 

ife  (2k+l) S 


Therefore  [1  -  2  S]  t(s)  =  y 

k=0 


1 

(2k+l)s 


for  Re(s)  ^  1  +  <$ ,  where  6  is  positive.  We  may  calculate  c(3/2)  and 
thus  obtain  the  sum  to  our  series.  The  integral  of  (B.48)  is  performed 
numerically  by  the  trapezoidal  rule. 
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OTHER  METHODS 


The  following  brief  account  will  be  outlined  with  respect  to  some  of  the 
above  mentioned  solution  procedures. 


Initially  Yuan8-? 


employed  the  simplified  Donnell  equation  of  the  form 
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to  obtain  a  solution  of  the  infinitely  long  cylinder  subject  to  two  equal  and 
opposite  point  loads.  The  loading  function  q(x,s)  was  resolved  in  a  product  of 
a  Fourier  integral  axially  and  a  cosine  Fourier  series  peripherally.  The 
loading  function  q(x,s)  in  radial  load  is  even  in  both  x  and  s  variables.  In 
their  follow-up  paper  Ting  4  Yuan8-®  used  the  complete  Donnell  equation,  which 
can  be  written  as  follows 
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They  dealt  with  boundary  conditions  other  than  simply  supported.  They  employed 
the  infinite  length  cylinder  solution  as  a  particular  integral.  Boundary 
conditions  at  the  ends  were  satisfied  by  obtaining  proper  solutions  to  the 
homogeneous  differential  equation. 

Finally  Yuan  &  TingB-9  employed  the  more  accurate  Fluegge  equation 
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(B.  53) 


The  method  of  images  was  employed  to  obtain  the  radial  displacement  of  a  simply 
supported  cylinder  under  the  applied  load  as 


wEh  =  6(l-v2)  [Rl2  |  1  COSh(Bl  I)'  cos  (al  |) 
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where  all  parameters  are  defined  in  Reference  B-9  or  in  main  text. 


(B. 54) 
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APPENDIX  C 

SYMMETRY  OF  VLASOV  SHELL  THEORY 


The  following  statements  are  made  to  establish  an  important  point  as  far  as 
the  form  of  the  used  shell  equations  is  concerned.  It  concerns  the  symmetry  of 
the  coefficient  matrix. 

Written  in  matrix  notation,  but  with  their  respective  sign  convention 
maintained,  the  thin  shell  Timoshenko  equations^-*-  used  in  Biljaard's 
solution, ^  ^  for  radial  load  pr  acquire  the  form  used  in  Equation  (C.l). 

Notice  the  asymmetry  of  the  coefficient  matrix  in  terms  C3  2  and  ^2  3* 

Fluegge's  equationsc-3  for  the  same  problem  in  his  notation  and  sign 
convention  are  shown  in  Equations  (C.2)  and  (C.3).  The  coefficient  matrix  is 
fully  symmetrical. 

The  same  is  true  with  Vlasov  equations .C-4,C-5  They  are  symmetrical  and 
can  be  displayed  as  shown  in  Equation  (C.4).  The  important  point,  which  more 
often  than  not  is  missed,  is  that  the  reciprocity  law  of  Betti  does  not  hold  for 
Timoshenko  equations . c_^  Within  the  limits  of  shell  theory,  however,  it  does 
not  restrict  their  application  in  the  "engineering"  sense.  Finally,  it  should 
be  clear  that  Fluegge,C-d  Vlasov ,C-^  or  other  types  of  shell  equations 
ought  to  be  preferred  due  to  symmetry  considerations. 
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